EXPLOSION  GENERATED  SEISMIC  WAVES 


M f* 
•^1 


M.  Nafi  Toksdz 

Department  of  Earth  and  Planetary  Sciences 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 


CONTRACT  NO.  F19628-68-C-0043 
Project  No.  7639 
Task  No.  763908 

Work  Unit  No.  76390801 


FINAL  REPORT 

Period  Covered:  1  October  1967  -  30  Septmeber  1971 


31  December  1971 


Contract  Monitor:  Ker  C.  Thomson 
Terrestrial  Sciences  Laboratory 


D  C 

lEGgrpnn 

y  KAY  a  W?  ,jjj 

JteEFU  biyJ 

C 


Approved  for  public  release;  distribution  unlimited. 


Prepared  for 


AIR  FORCE  CAMBRIDGE, RESEARCH  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
BEDFORD,  MASSACHUSETTS  01730 

NAT  ION/' I  irr.MNICAl 
INf  OH  MAT  ION  SfPVICf 


!ft,  asaect^i^  n  awi^M  H^ajjkiiM'UJii  [■jiUMOr^irii.  j^im: 


UNCLASSIFIED 


Sjruntjr^jl***»fic*tior^ 


-  139  - 


DOCUMENT  CONTROL  DATA  •  R  &  D 

fSacunty  clattillcalion  ot  title.  body  of  abatr*.  •  mnd  indexing  mnnotntson  iiioa#  to#  entered  »nw  the  overall  report  la  c tmaallledj 


ORiGin*  TING  activity  (Corporate  ruthor) 

Massachsuetts  Institute  of  Technology 
Department  of  Earth  and  Planetary  Sciences 
Cambridge,  Massachusetts  02139 


7j».  Btcont  sccjri’v  c l Ainr ic * tion 

UNCLASSIFIED 


2b.  CROUP 


3  REPORT  TITLE 

EXPLOSION  GENERATED  SEISMIC  WAVES 


4  descriptive  NOTES  (Ty*  e  ot  report  mnd  Incluatve  dmtma) 

Scientific. _ Final:  1  October  1967-30  September  1971 

*  ao  THORlSt  ffirit  name,  middle  initial.  Imat  nmma) 

M.  Nafi  Toksflz 


Approved  -• 
02-18-72 


t  REPORT  DATE 

31  December  1971 

•*.  COM  TRAC  T  OR  GRANT  NO 

F19628-68-C-0043 

6.  ahojcct  no  f  Task  Work  Unit  No. 

7639-08-01 

c 

DoD  Element  62101F 
_ .LQD— SubfileuLfcnt _ 681M.Q _ 

j  1C  DISTRIBUTION  STATEMENT 

A  -  Approved  for  public  release; 


la.  TOTAL  NO  OP  PAGES 


141 


lb.  NO  OF  RErj 


75 


9m.  ORIGINATOR'S  REPORT  NUUBCRISI 


9b.  otmCR  REPORT  NOISI  (Any  othmr  numbmra  that  may  bm  aealgned 
thia  report) 

AFCRL-72-0012 


distribution  unlimited. 


1.  JUPPLCUtNUSV  NOTE* 

TECH,  OTHER 


1J  ABSTRACT 


13  SPONSORING  MILItAR*  ACTIVITV 

Air  Force  Cambridge  Research  Labora 
tories, (LW) 

L.  G.  Hanscom  Field 
Bedford.  Massachusetts  01730 


The  mechanisms  of  generation  of  SH  and  Love  waves  by  under¬ 
ground  nuclear  explosions  are  investigated  theoretically  and  experi¬ 
mentally.  An  explosion  in  pte-s tressed  medium  is  formulated  as  a 
composite  source  consisting  of  an  isotropic  explosion  and  a  double¬ 
couple  representing  the  tectonic  strain  release.  From  the  radiation 
patterns  of  Rayleigh  waves  and  Love-to-Rayleigh  wave  amplitude  ratios 
relative  strength  of  the  strain  release  component  is  determined.  It 
is  found  that  relative  amplitudes  of  Love  waves  and  tectonic  strain 
energy  release  strongly  depend  on  the  medium  properties  in  the  immedi¬ 
ate  vicinity  of  the  explosion.  For  "harder"  media  such  as  granite, 
energy  release  is  much  greater  compared  to  "softer"  media  such  as  al¬ 
luvium  or  salt.  The  strain  energy  release,  as  best  as  can  be  deter¬ 
mined  from  laboratory  experiments,  is  primarily  due  to  the  relaxation 
of  the  pre-stress  field  around  explosion-generated  and  extended  cracks 
and  faults. 


I 

t 


DD  :r„1473 


UNCLASSIFIED 

Srrufity 


UNCLASSIFIED 


Security  Classification 


-  140  - 


Underground  explosions 
Source  mechanisms 
Tectonic  strain  release 
Radiation  pattern 


EXPLOSION  GENERATED  SEISMIC  WAVES 

M.  Nafi  Toksflz 

Department  of  Earth  and  Planetary  Sciences 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 


CONTRACT  NO.  F19628-68-C-0043 
Project  No.  7639 
Task  No.  763908 

Work  Unit  No.  76390801 


FINAL  REPORT 


Period  Covered: 


1  October  1967  - 

31  December  1971 


1971 


Contract  Monitor:  Ker  C.  Thomson 
Terrestrial  Sciences  Laboratory 


Approved  fcr  public  release;  distribution  unlimited. 


Prepared  for 

AIR  FORCE  CAM3RIDGE, RESEARCH  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 


BEDFORD,  MASSACHUSETTS  01730 


-  11  - 

ABSTRACT 

The  mechanisms  of  generation  of  SH  and  Love  waves  by 
underground  nuclear  explosions  are  investigated  theoretically 
and  experimentally.  An  explosion  in  pre-stressed  medium  is 
formulated  as  a  composite  source  consisting  of  an  isotropic 
explosion  and  a  double-couple  representing  the  tectonic  strain 
release.  From  the  radiation  patterns  of  Rayleigh  wave;,  and 
Love-to-Rayleigh  wave  amplitude  ratios  relative  strenath  of 
the  strain  release  component  is  determined.  It  is  found  that 
relative  amplitudes  of  Love  waves  and  tectonic  strain  energy 
release  strongly  depend  on  the  medium  properties  in  the  im¬ 
mediate  vicinity  of  the  explosion.  For  "harder"  media  such 
as  granite,  energy  release  is  much  greater  compared  to  "softer" 
media  such  as  alluvium  or  salt.  The  strain  energy  release, 
as  best  as  can  be  determined  from  laboratory  experiments,  is 
primarily  due  to  the  relaxation  of  the  pre-stress  field 
around  explosion-generated  and  extended  cracks  and  faults. 
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INTRODUCTION 

The  source  mecnarism  of  underground  nuclear  explosions 
lias  been  of  great  interest  to  seismologists  since  it  was  first 
observed  that  a  pure  explosive  source  alone  could  not  adequ¬ 
ately  account  for  the  observed  seismic  radiation.  The  exis¬ 
tence  of  prominent  SH  and  Love  waves  in  the  ~eismic  records 
from  these  events  was  a  prime  reason  for  this  conclusion  and 
indicated  that  the  explosions  were  probably  causing  the  release 
of  some  tectonic  strain  (Press  and  Archambeau,  1962).  Early 
studies  by  Brune  and  Pomeroy  (1963),  Aki  (1964),  and  Toksflz 
et  al  (1963)  found  that  a  double-couple  type  source  combined 
with  a  pure  explosion  could  account  for  the  observed  seismic 
radiation  from  several  events. 

In  this  report  we  summarize  the  work  we  have  done  on  a 
number  of  nuclear  explosions  at  the  Nevada  Test  Site  and  in 
other  areas.  This  report  will  consist  of  three  parts,  each 
of  which  is  essentially  self-contained.  There  is,  therefore, 
a  certain  amount  of  repetition,  such  as  in  the  theory,  some 
of  the  tables,  and  in  the  references  contained  in  the  report. 
Part  I  deals  with  the  source  mechamisms  determined  for  a  large 
number  of  underground  explosions  and  the  comparison  of  these 
results  to  other  data  connected  v/ith  these  events.  Part  II 
consists  of  a  detailed  analysis  of  one  nuclear  event  (Silby) 


and  •  report  concerning  laboratory  experiments  on  explosive 
sources.  The  combination  of  these  wc  types  of  analyses 
helps  to  clarify  the  problem  of  strain  energy  release  in  pre- 
st.ressed  media.  Part  III  is  concerned  with  two  nuclear 
events  (Pile  Driver  and  Greeley)  which  released  a  significantly 
ru  jher  proportion  of  tectonic  strain  energy  than  other 
v-x,. )  .  ions. 
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Introduction 

The  source  mechanisms  have  been  determined  for  a  number 
of  underground  explosions  at  the  Nevada  Test  Site,  in  the 
Soviet  Union,  and  in  other  .{testing  sites.  We  have  followed 
the  method  of  Tokstiz  et  al.  (1965)  of  considering  the  source 
to  be  a  superposition  of  a  symmetrical  explosion  and  a  hori¬ 
zontal  double-couple.  The  Love  wave  to  vertical  component 
Rayleigh  wave  amplitude  ratio  from  such  a  source  is  then  given 


F  kf  AL  Cos-ae 

( t  +  f  Sin  &e)  ^  CU'M~) 


where  assuming  the  source  time  functions  for  both  the  explo¬ 
sions  and  double-couple  and  the  attenuation  coefficients  for 
Love  and  Rayleigh  waves  are  the  same,  kL  and  kR  are  the 
Love  and  Rayleigh  wave  numbers ,  Aj  and  AR  are  the  medium 
responses  for  Love  and  Rayleigh  waves,  uq  and  wq  are  the 
components  of  particle  velocity  at  the  surface,  F  is  the 
relative  strength  of  the  doub-e-eouple ,  and  9  is  azimuth 
of  the  fault  plane.  We  determined  the  two  parameters  F  and 
6  for  each  explosion  by  fitting  the  observed  surface  wave 
data  to  the  above  formulation. 
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Data 

Tables  la  and  lb  list  the  nuclear  explosions  we  have 
studied.  Included  are  seven  events  previously  analyzed  by 
Tokstiz  et  ax,  (1965),  Toksfiz  and  Kehrer  (1971),  and  Toksfiz 
(1967) .  The  surface  wave  data  was  obtained  from  stations 
throughout  North  America  i.n  the  case  of  the  American  events. 
These  included  stations  of  the  World  Wide  Standard  Network 
(WWSS) ,  the  Long  Range  Measurements  Network  (LRSM) ,  the  Canadian 
Standard  Seismograph  Network  (CSS) ,  and  two  additional  stations 
at  Berkeley  and  Pasadena  which  provided  long-period  data.  Figure 
1  shows  the  distribution  of  recording  stations  in  North  America. 
However,  not  all  stations  provided  useful  data  for  each  event. 
For  the  presumed  Soviet  events,  data  was  obtained  from  the  WWSS 
stations  in  Europe  and  Asia. 

The  film  records  of  those  events,  for  which  the  surface 
wave  data  was  of  sufficient  guality,  were  manually  digitized. 

The  resulting  data  was  then  band-pass  filtered  to  reduce  noise 
and  the  two  horizontal  components  were  rotated  to  radial  ar.d 
tangential  directions  with  respect  to  the  epicenter.  In  this 
way  essentially  pure  Love  waves  were  obtained  on  the  tangential 
component.  Figure  2  is  an  example  of  this  analysis  for  the 
explosion  Boxcar  at  Mould  Bay.  Canada.  The  three  lower  traces 
are  plotted  from  digital  data  and  show  the  unfiltered  north- 
south,  east-west,  and  vertical  components.  The  three  upper 


traces  show  the  filtered  data,  where  the  horizontal  components 
have  been  rotated  to  radial  and  tangential  directions. 

The  tangential  and  the  vertical  components  were  then 
Fourier-analyzed  to  yield  the  amplitude  spectra.  The  proper 
time  windows  for  this  analysis  were  determined  by  examining 
the  filtered  trace  at  each  station  individually  and  noting 
the  onset  and  duration  of  the  wave  trains. 

Love  wave  to  Rayleigh  wave  amplitude  ratios  were  deter¬ 
mined  by  computing  the  spectral  ratio  of  the  tangential  com¬ 
ponent  to  the  vertical  component  over  the  period  range  8  to 
100  seconds.  In  general,  the  spectral  ratios  were  somewhat 
oscillatory  and  thus  the  ratio  at  any  one  particular  period 
could  not  be  expected  to  be  reliable  at  all  stations.  The 
period  range  was,  therefore,  divided  into  three  bands,  from 
9  to  1j  seconds,  from  15  to  22  seconds,  and  from  22 
t.o  30  seconds.  The  average  value  of  the  ratio  in  each  band 
wv.s  then  taken  to  represent  the  ratio  ov^r  that  banu. 

An  automatic  error  scheme  (Kehrer,  1969)  was  used  to  deter¬ 
mine  the  best  fit  of  Equation  [1]  to  the  observed  data.  In 

■»  t  |  y<y 

applying  Equation  [1],  the  term  AR(~)  was  taken 

o 

to  be  a  constant  at  all  stations.  7m  error  term  E  was  then 
formed  between  the  observed  and  the  theoretical  ratios  for 
combinations  of  F  and  0. 


i  f  /N 

1-vFj  W  ' 
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The  jth  combination  of  F  and  6  which  fits  the  data  best 
will  be  that  which  minimizes  E.  L^/ft ^  is  the  experimentally 
measured  Love  wave  to  Rayleigh  wave  amplitude  ratio  at  a  par 
ticuiar  station  i,  S  is  a  constant,  and  N  is  the  number 
of  stations.  The  values  of  E  were  contoured  on  a  Stromberg- 
Carlson  4020  grind  for  the  three  period-bands  of  each  event. 
Figure  3  shows  these  results  for  the  explosion  Faultless.  The 
minima  indicate  the  orientation  of  the  best  fitting  right- 
lateral,  vertical,  strike-slip  fault  as  measured  clockwise 
from  the  north  at  the  source.  A  comparison  of  the  error  plots 
over  the  three  bands  in  Figure  3  reveals  the  frequency  depen¬ 
dence  of  r.  The  orientation  8  of  the  double-couple  remains 
essentially  unchanged  over  the  three  bands  while  F  is  greatest 
at  long  periods  and  smallest  at  short  periods.  This  relation 
ship  was  observed  for  nearly  all  events  and  indicates  that  the 
source  of  Love  waves  is  less  efficient  at  short  periods  than 
the  Rayleigh  wave  source.  A  similar  phenomenon  was  noted  oy 
Akj.  et  al.  (1969)  for  body  waves  frcm  the  Benharn  event.  Figure 
4  shows  the  relationship  between  F  and  the  period  for  several 
events. 

Two  other  matrices,  in  addition  to  the  one  described,  were 
investigated  to  test  the  best  fit  of  the  data  to  theoretical 
source  configurations.  The  one  consisted  of  normalizinq  the 
difference  between  the  observed  and  theoretical  ratios  by  tne 
observed  ratio  at  each  station.  This  reduced  the  weights  of 
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large  observed  values  on  the  error  term#  with  the  effect  that 
the  F  value  was  reduced.  The  second  method  consisted  or 
normalizing  by  the  theoretical  ratio.  This  had  the  effect 
of  increasing  the  F  value.  In  general,  however,  the 
changes  in  F  were  small  and  the  orientation  remained  essen¬ 
tially  unchanged.  The  values  of  F  reported  here  are  those 
determined  by  the  non-normalized  matrix  with  the  exception  cf 
Creeley.  In  the  case  of  Greeley,  the  absolute  minimum  in  the 
error  plot  for  the  15  to  22  second  period  range  occurred 
for  F  =  0.9  and  0  =  175°.  However,  this  configuration 
could  not  account  for  the  consistently  large  values  of  the 
ratio  due  north  of  the  NTS.  Over  the  22  to  30  second 
range,  the  indicated  miminum  occurred  at  the  much  larger  F 
value  of  1.6.  Indeed  this  value  appeared  to  fit  the  data 
best  over  the  15  to  22  second  range  as  well  and  it  was 
therefore  taken  to  represent  the  strength  of  the  double-couple 
for  Greeley.  Furthermore,  an  investigation  of  Rayleigh  wave 
polarity  (Toksttz  and  Kohrer,  1971)  indicated  that  the  F  value 
for  Greeley  was  greater  than  one. 

A  summary  of  the  best  fitting  F  and  9  in  the  15  to 
22  second  period  range  for  each  event  is  given  in  Tables  la 
and  lb.  For  the  explosions  of  previous  studies  and  for  Buff, 
Benham,  and  the  presumed  Soviet  events  of  3-20-66,  3-26-67, 
and  10-21-67  for  which  the  surface  wave  data  was  not  of  suf¬ 
ficient  quality  to  permit  the  calculation  of  spectral  ratios, 


the  parameters  F  and  o  were  determined  from  mea¬ 
surements  of  the  maximum  amplitudes  of  Lo-'e  and  Rayleigh 
waves.  Included  in  Table  I  is  the  ratio  of  the  surface  #ave 
energy  of  the  double-couple  to  the  energy  of  the  explosion. 
This  ratio  (Toksdz  et  al,  1965)  is  closely  approximated  by 

Etect/Exp  -  </3  F?'- 

Several  conclusions  can  be  drawn  from  Table  I.  The 
close  agreement  in  fault  orientation  of  most  of  the  events  at 
the  Nevada  Test  Site  indicates  that  a  regional  stress  field 
is  probably  controlling  the  strain  release.  The  large  vari- 
Lion  v’ithin  the  two  Soviet  testing  areas  may  be  due  to  a  lack 
of  data  in  the  case  of  three  of  the  events.  The  dependence 
of  F  and  hence  the  amount  of  stress  release,  on  lithology 
was  previously  noted  by  Toksflz  (1967)  and  by  Toksfiz  and  Kehrer 
(1971) .  Table  I  provides  additional  evidence  that  the  multi¬ 
polar  source  component  is  larger  in  more  competent  rock.  Thus 
strain  release  from  a  particular  explosion  is  dependent  on 
the  rock  medium  and  the  regional  state  of  strain.  The  rela¬ 
tively  low  F  values  for  all  of  the  presumed  Soviet  events 
are  consistent  with  what  is  known  about  the  two  Soviet  testing 
areas. 

The  north-northwest  orientation  of  the  fault  planes  of 
most  of  the  Nevada  Test  Site  explosions  is  in  good  agreement 
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with  lineaments  and  dominant  fault  patterns  in  the  area. 

Figure  5  shows  the  Yucca  Fait  portion  of  the  Nevada  Test  Site 
with  the  source  mechanisms  of  the  explosions  studied  there. 
Included  are  the  major  natural  faults  and  some  of  the  fractures 
produced  by  the  events.  The  chief  feature  of  Yucca  Flat  is 
the  Yucca  Fault  which  trends  northward  through  alluvium  and  is 
of  recent  age.  Although  most  observed  displacements  across 
the  Yucca  Fault  are  vertical  following  nearby  explosions, 
numerous  en  echelon  fractures  along  the  fault  were  produced 
in  the  alluvium  by  Corduroy  (Barosh,  1968).  Such  an  en  echelon 
pattern  is  generally  related  to  strike-slip  movement  in  the 
underlying  basement.  The  trend  of  the  pattern  in  this  case 
suggests  right-lateral  slippage. 

Figure  6  shows  the  Pahute  Mesa  portion  of  the  Nevada 
Test  Site  with  the  faults,  explosion-produced  fractures,  and 
the  events  studied.  The  main  structural  feature  of  PahvLe 
Mesa  is  the  Silent  Canyon  Caldera,  which  encloses  the  five 
events.  Many  normal  faults,  strikinq  north-northwest,  cut 
the  thick  sequence  of  Tertiary  volcanic  rock.  Recent  natural 
movement  along  some  of  these  faults  has  been  inferred  (McKeown 
et  al.,  1966).  For  Benham  and  Boxcar,  studies  of  aftershock 
distribution  (Hamilton  and  Healy,  1969;  Ryall  and  Savage,  196S) 
as  well  as  studies  of  faulting  and  fracturing  caused  by  the 
explosions  (Buckman,  1969;  McKeown  and  Dickey,  1969)  have  beer, 
published.  Figure  7  shows  the  relationships  of  these  two  sets 
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of  data  with  the  source  mechanisms  determined  for  the  two 
events.  In  the  case  of  Benham,  the  correlation  between  the 
orientation  of  the  double-couple  and  the  trend  of  the  after¬ 
shock  pattern  as  well  as  the  orientation  of  faults  and  frac¬ 
tures  is  quite  good.  Both  vertical  and  right-lateral  move¬ 
ments  on  the  nearby  Boxcar  fault  were  initiated  by  Benham. 

Hamilton  and  Healy  (1969)  found  that  the  aftershocks  of 
Benham  were  concentrated  in  two  patterns.  Those  to  the  north¬ 
west  of  ground  zero  follow  a  northeasterly  trend  while  those 
to  the  west  and  southwest  trend  north-south.  Fault  plane 
solutions  from  P  waves  of  events  in  the  first  group  indicate 
dip-slip  fault  movement.  Solutions  for  the  second  group  in¬ 
dicate  right-lateral  strike-slip  motion  in  a  northerly  direction. 
These  two  patterns  were  further  confirmed  in  a  study  of  smaller 
Benham  aftershocks  by  Stauder  (1971).  The  solutions  for  the 
second  group  of  aftershocks  are  similar  to  the  orientation  of 
the  double-couple  component  for  Benham  (and  other  events)  as 
determined  in  this  study.  The  Benham  and  Boxcar  aftershocks 
when  considered  together  (see  Figure  7)  indicate  the  continuity 
of  the  distribution  pattern.  Northeast-trending  Benham  after¬ 
shocks  overlap  with  those  of  Boxcar. 
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Figure  Captions 

Figure  1.  Stations  recording  surface  waves  from  explosions 
in  the  continental  United  States. 

Figure  2.  Filtered  (upper  three  traces)  and  unfiltered 

(lower  three  traces)  long-period  components  from 
the  Boxcar  explosions  at  Mould  Bay,  Canada  (MBC) 
plotted  from  digital  data  at  the  same  gain.  The 
two  filtered  horizontal  components  have  been 
rotated  to  radial  and  tangential  directions  to 
separate  Love  and  Rayleigh  waves. 

Figure  3.  Contour  plots  of  the  deviations  of  combinations 
of  part  double-couple  and  fault  plane  azimuth 
from  experimental  Love  to  Rayleigh  wave  ratios 
for  three  period  ranges:  a)  9  to  15  sec, 
b)  15  to  22  sec,  and  c)  22  co  jC  sec, 
for  the  Faultless  event. 

Figure  4.  The  period  dependence  of  F,  the  past  double¬ 
couple,  for  five  nuclear  expljsions. 

Figure  5.  Map  of  Yucca  Flat  showing  the  fault  plane  orienta¬ 
tion  of  the  explosions  studied,  explosion -produced 
fractures,  and  major  natural  faults  in  the  area 
(geology  after  Fernaid  et  al,  1968). 
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Figure  6.  Map  of  Pahute  Mesa  showing  the  fault  plane 

orientations  of  the  explosions  studied,  explosion- 
produced  fractures,  and  natural  faults  in  the 
area  (faults  after  Orkild  ot  al,  1969). 

Figure  7.  Source  mechanisms  determined  for  the  Boxcar  and 
Benham  events  compared  to  natural  faults  and 
explosion-produced  fracturing,  and  aftershock 
distribution.  Aftershock  distribution  for  Benham 
(Hamilton  and  Healy,  1969)  is  that  of  the  larger 
aftershocks  (magnitude  3  to  4.2).  Aftershocks 
from  Boxcar  (Ryall  and  Savag.-,  1969)  are  those 
occurring  in  the  first  two  days  after  the  event. 
The  circles  are  of  radius  5  km  centered  on  the 
ground  zero  of  the  explosion. 
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The  mechanisms  of  generation  of  seismic  waves  by  an 
explosion  in  prestressed  media  are  studied  using  both  field 
seismograms  and  controlled  laboratory  experiments.  LRSM 
seismograms  from  the  underground  nuclear  explosion  Bilby  are 
analyzed  to  determine  the  source  parameters  from  the  radia¬ 
ted  Love  and  Rayleigh  waves.  From  the  normalized  amplitudes 
of  Rayleigh  waves  as  well  as  the  Love-Ray leigh  amplitude  ratios, 
a  composite  source  consisting  of  an  isotropic  explosion  and 
a  double-couple  is  synthesized  for  the  explosion  and  the 
associated  tectonic  strain  release.  From  Bilby  and  ouher 
explosions  studied  by  similar  techniques,  it  is  found  that 
the  tectonic  strain  energy  release  strongly  depends  on  the 
medium  properties  in  the  immediate  vicinity  of  the  explosion. 

For  "harc>r"  media  (such  as  granite)  the  tectonic  strain 
energy  release  and  the  relative  amplitude  of  Love  waves  are 
significantly  higher  than  fcr  softer  media  such  as  alluvium. 
Source  time  functions  of  Love  waves  associated  with  the 
explosions  are  closer  to  time  functions  of  earthquakes  than 
to  those  of  explosions. 

The  mechanisms  of  the  pre-existing  strain  energy  release 
by  explosive  sources  are  studied  in  two  separata  laboratory 
experiments.  In  a  one-dimensional  experiment  where  an  ex¬ 
plosive  source  ii.  detonated  in  a  rod  stressed  in  torsion,  the 
S-vave  amplitudes  are  found  to  be  linearly  proportional  to 
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pre-strain.  In  the  second  experiment,  radiation  of  seismic 
waves,  and  the  near  source  phenomena  of  explosive  sources 
in  pre-stressed  plates  are  studied  by  photoelastic  as  well 
as  strain  gauge  observations.  The  generation  of  S-waves 
is  greatly  enhanced  by  the  pre-stress  condition.  It  is 
found  that  extended  cracking  (faulting)  occurs  along  direc¬ 
tions  determined  by  the  prestress  field.  The  transverse 
(SH)  waves  are  generated  primarily  by  the  relaxation  of  the 
stress  field  along  these  cracks.  The  explosion-generated 
cavity  alone  could  not  account  for  the  radiated  transverse 
seismic  energy. 
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I .  INTRODUCTION 

Seismic  surface  and  body  waves  generated  by  underground 
explosions  and  recorded  at  distant  stations  quite  often  show 
characteristics  differing  from  those  of  an  ideal  explosive 
point  source  in  a  homogeneous  medium.  Most  notable  among 
these  differences  are  the  presence  of  Love  waves  and  the 
azimuthal  asymmetry  of  the  Rayleigh  wave  radiation  patterns. 
Through  numerous  studies  it  has  been  shown  that  the  Love  waves 
are  produced  at  the  source  region,  although  the  exact  mecha¬ 
nism  of  their  generation  still  cannot  be  clearly  defined. 

The  purpose  of  this  paper  is  to  present  and  interpret  data, 
obtained  from  underground  nuclear  explosions  and  controlled 
laboratory  experiments,  on  seismic  waves  generated  by  explo¬ 
sive  sources  in  stressed  media. 

There  have  been  a  number  of  studies  dealing  with  the 
radiation  patterns  and  the  generation  of  SH-type  seismic  waves 
by  underground  nuclear  explosions  (Press  and  Archambeau,  1962; 
Brune  and  Pomeroy,  1963;  Aki,  1964;  Toksttz  et  al . ,  1964; 

Toksdz  et  al.,  1965;  Toksfiz,  1967;  Kehrer,  1969;  Molnnr  ct 
.  al. ,  1969;  Tsai  and  Aki,  1970).  In  addi.tion,  there  arc  a 
number  of  previous  theoretical  studies  on  the  problems  o£ 
seismic  wave  generation  by  explosions  (Zvolinskii,  1950;  Bishop, 
1963;  Alverson,  1964;  Cisternas,  1964;  Butkovich,  1965; 
Archambeau,  1968).  Difficulties  in  describing  wave  propagation 
from  nuclear  explosions  arise  because  of  very 
high  stresses  involved.  In  the  i:v  lia'.  •  'J  ' 
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extremely  important  to  the  understanding  of  seismic  wave 
generation,  yet  this  region  is  extremely  difficult  to  de¬ 
scribe  mathematically. 

In  the  description  of  this  complex  problem,  difficulties 
arise  both  in  mathematical  formulation  and  from  inadequate 
knowledge  of  the  behavior  of  geologic  materials  under  shock 
loading.  In  the  zone  immediately  around  the  explosion,  a 
relatively  small  volume,  hydrodynamic  compressible  flow 
equations  can  be  used  without  major  assumptions  or  approxima¬ 
tions.  In  the  zones  where  the  rocks  are  crushed  and  cracked, 
however,  neither  fluid  nor  elastic  behavior  can  be  assumed, 
and  a  semiempirical  treatment  based  on  approximations  and 
available  data  must  be  followed.  If  the  explosion  is  placed 
in  a  medium  where  tectonic  stresses  exist,  other  complica¬ 
tions  arise  from  the  stress  relaxation  around  the  cavity,  the 
crushed  zone,  and  the  zone  of  numerous  shock  induced  cracks. 

Controlled  laboratory  experiments  can  provide  some  of 
the  information  needed  for  the  understanding  of  seismic  wave 
generation  by  explosions  in  pre-stressed  media.  Some  experi¬ 
ments  have  been  conducted  in  the  field  using  explosive  sources 


in  soil  (Kisslinger  et  al . ,  1961).  Other  model  experiments  have 
been  carried  out  in  the  laboratory  using  explosive  sources  in 
pre-stressed  plates.  In  these  tests,  the  radiat’.on  patterns  of 
seismic  compressional  and  shear  waves  were  observed  with 
various  transducers  and  oscilloscopes  (Kim  and  Kisslinger,  1967; 


Kisslinger  and  Gupta,  1963).  Most  recently,  photoelastic 
techniques  have  been  used  to  observe  wave  generation  and  pro¬ 
pagation  in  pre-stressed  media  (Thomson  et  al . ,  1969). 

In  this  paper  we  will  treat,  in  three  steps,  the  problem 
of  mechanisms  of  seismic  wave  generation  by  explosive  sources 
in  pre-stressed  media.  First,  we  describe  the  radiation 
pattern  of  surface  waves  from  a  typical  underground  explosion 
at  the  Nevada  Test  Site.  We  determine  its  source  mechanism 
by  using  an  amplitude  equalization  scheme.  In  part  two, 
we  describe  the  results  of  laboratory  experiments  in  which 
explosive  sources  were  detonated  in  plates  ur».»er  different 
pre-stress  conditions.  Here  the  mechanisms  of  seismic  wave 
generation,  crack  formation,  and  source  complications  were 
observed  using  dynamic  photoelasticity  and  high  speed  photo¬ 
graphic  techniques.  Finally,  we  interpret  the  field  obser¬ 
vations  in  light  of  laboratory  findings. 

II.  FIELD  OBSERVATIONS 

Seismic  waves  from  a  large  number  of  the  U.  S.  underground 
explosions  have  been  well-recorded  by  LRSM  (Long  Range  Seismic 
Measurements) ,  WWSS  (World-Wide  Seismic  System)  and  other 
seismograph  stations  in  North  America.  LRSM  stations  provide 
the  best  recordings  of  surface  waves  in  the  period  range  of 
10  to  50  seconds.  Most  o£  the  large  explosions  detonated  in 
relatively  hard  geologic  media  (granite,  tuff)  generate  Love 
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waves  in  addition  to  Rayleigh  waves.  Two  examples  of  Love 
waves  from  explosions  are  shown  in  Figures  1  and  2 .  Explo¬ 
sions  in  loose  alluvium  and  salt  domes,  and  collapse  events 
following  explosions,  do  not  generate  a  significant  amount 
of  Love  waves.  The  main  body  of  this  study  is  devoted  to 
the  understanding  of  the  generation  of  these  transverse 
waves.  Some  sources  that  must  be  considered  are: 

(1)  relaxation  of  the  pre-stressed  medium  around  the  explo¬ 
sion  generated  cavity,  (2)  triggering  of  an  earthquake  by 
the  explosion,  (3)  radiation  from  explosion  induced  cracks 
in  the  medium,  and  (4)  a  combination  of  all  of  these. 

In  order  to  resolve  this  problem  it  is  useful  to  first 
determine  the  radiation  patterns  of  Rayleigh  and  Love  waves. 

A  group  of  explosions  have  been  studied  using  the  amplitude 
equalisation  method  (Toksflz,  e_t  aT .  ,  1964,  1965;  Toksdz,  1967; 
Kehrer,  1969) .  Here  we  illustrate  the  method  with  data  from 
the  Bilby  explosion. 
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II. 1.  Radiation  of  Seismic  Waves  from  the  Bilby  Explosion. 

The  nuclear  explosion  Bilby  was  detonated  at  the  Nevada 
Test  Site  on  13  September  1963  at  17:00:00  GMT.  The  shot  was 
placed  at  a  depth  of  700  m.  in  tuff.  The  computed  equivalent 
magnitude  was  m^  =  5.8.  A  collapse  event  (m^  =  4.5)  followed 
the  explosion  at  17:31:20.5  GMT. 

Bilby  generated  Rayleigh  and  Love  waves  which  were  well 
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recorded  at  a  number  of  LRSM  stations.  The  station  distri¬ 
bution  is  shown  in  Figure  3.  The  Love  waves  could  be  iden¬ 
tified  easily  since  the  long-period  hori2ontal  instruments 
were  generally  oriented  in  radial  and  transverse  directions 
relative  to  the  source.  Furthermore,  fora  continental  path. 
Love  waves  have  higher  group  velocities  than  Rayleigh  waves 
in  the  period  range  of  interest. 

Radiation  Pattern  of  Rayleigh  Waves .  Under  ideal  conditions 
the  radiation  pattern  can  be  >btained  by  correcting  the 
observed  amplitudes  for  instrument  response,  geometric  spread¬ 
ing  and  attenuation  effects.  In  practice,  however,  the  geo¬ 
logic  conditions  are  complicated  and  the  absolute  amplitude 
method  is  not  very  effective  (Toksfiz  and  Clermont,  1967). 

Thus  a  method  of  normalization  must  be  considered.  The 
source  and  the  propagation  factors  cannot  be  isolated  from 
each  other.  Therefore,  we  must  consider  methods  in  which 
one  of  the  variables  is  held  constant  while  the  effect  of 
the  other  is  investigated.  We  used  two  methods  for  the 
study  of  relative  amplitudes  to  determine  the  radiation  pat¬ 
tern  from  Bilby.  In  the  first,  we  normalized  the  amplitudes 
of  the  Rayleigh  waves  generated  by  the  explosion  to  those 
of  the  collapse  event  that  followed.  In  the  second  we  used 
the  ratio  of  the  Love  and  Rayleigh  wave  amplitudes  to  obtain 
a  source  function. 

The  collapse  of  the  cavity  (formed  by  the  explosion)  seems 
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to  provide  an  excellent  reference  for  normalization.  Whan 


the  amplitude  ratios  of  the  cxplosion-and  collapse—  generated 


waves  are  taken ,  the  propagation  and  instrument  effects  com¬ 


pletely  cancel  out,  and  the  ratio  directly  reflects  the  source 


effects.  In  previous  studies,  it  was  found  that  the  radia¬ 


tion  pattern  from  the  post-explosion  collapse  was  in  ‘general  more 


symmetric  than  the  pattern  from  the  explosion.  The  azimuthally 


uniform  amplitudes  as  well  as  the  absence  of  prominent  Love 


waves  were  evidence  of  this  radial  symmetry  (Toksfiz,  1967). 


If  we  assume  that  the  Bilby  collapse  had  a  radially  uniform 


radiation  pattern,  then  the  explosion/collapse  ratio  would 


indicate  if  the  explosion  was  a  symmetric  source. 


The  ratios  of  peak  amplitudes  of  the  Rayleigh  waves  are 


shov;n  in  Figure  4  as  a  function  of  azimuth.  The  peak  ampli¬ 


tudes  were  directly  read  from  the  long-period  records  of  the 


LRSM  stations.  In  a  few  instances  the  sr.ectral  ratios  were 


computed,  and  these,  on  the  average,  showed  a  similar  radiation 


pattern  as  the  peak  amplitudes. 


The  azimuthal  cov<  rage  in  Figure  4  is  far  from  being 


complete,  but  in  the  north  and  northeast  directions,  where 


there  is  a  sufficient  number  of  observations,  the  deviation 


from  a  uniform  ratio  is  clear.  If  we  assume  that  the  radial 


non-uniformity  of  the  radiation  pattern  from  the  explosion  was 


due  to  some  form  of  tectonic  complication  (such  as  the  relaxa¬ 


tion  of  the  medium  due  to  the  cavity  or  induced  rupture) ,  we 


can  include  this  effect  by  superimposing  a  multipolar  term 


■ 
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on  the  explosive  source.  Both  seismic  model  experiments  and 


theoretical  studies  indicate  that  at  large  distances  from  the 
source,  a  double-couple  type  source  function  is  a  good  repre¬ 
sentation  for  tectonic  strain  release  (Honda,  1962;  Burridge 
and  Knopoff,  1964;  Archambeau,  1968).  Then  the  displacements 
observed  at  a  distant  station  can  be  written  as  the  vectorial 
sums  of  those  due  to  an  explosive  source  and  to  a  double¬ 
couple  . 

Using  the  notation  of  Toksfis,  et  al.  (1965),  we  can  write 
the  far-field  expressions  for  the  Rayleigh  wave  ground  dis¬ 
placements  from  a  near-surface  explosive  source. 

c  /  U*  \ 

We{o)  =  k^/2  (  J  AR(ca)T(a))exp(-YRr)exp(i(wt-kRr-«J>t+3Tr/4)3 


Gi  /u*V 

U  (u)  »  * - k^*  (  ™  ]  A  (io)T(u)axp(-Ypr.)exp[i  (u)t-kftr-4)  -3it/4)  ] 

®  (2Ttr)H  R  \W0/ 


V  (a)  »  0 


W  («) ,  U  (so) ,  V  («)  are  the  vertical,  radial  and  tangential 
components  of  the  displacement,  kR  is  the  wave  number,  r  is 
the  radial  distance,  yr  is  the  Rayleigh  wave  attenuation  coef¬ 
ficient.  Ar  is  the  medium  response  for  Rayleigh  waves  due  to 
a  vertical  force,  U  and  WQ  are  the  components  of  particle 
velocity  at  the  surface.  T(a}  and  $fc{u))  are  the  amplitude  and 


4>c  -  tt  (Brune  and  Pomeroy,  1963;  Smith,  1963;  Toksfiz  et  al. , 
1964),  With  the  above  formulations  and  assumptions,  the 
theoretical  Rayleigh  wave  displacement  ratios  for  a  given 
explosion-collapse  pair  can  be  written  as 


(Ur) 

explosion  „  c, (1  .  F  sin  2e) 

(up) 

collapse 


(6) 


In  Figure  4  the  theoretical  curve  is  computed  using  (6) . 
Choosing  C' ,  F,  and  the ‘orientation  of  the  double-couple  to 
fit  the  data  best,  we  find  C*  =  12,  F  =  0.47,  and  the  refer¬ 
ence  direction  for  the  double-couple  principal  plane  (i.e, 
"fault  plane")  8  =  340°.  With  these  data  the  fit  is  good. 
The  above  figures  mean  that  the  explosion-generated  surface 
waves  were  12  times  larger  than  those  of  the  collapse,  and 
that  the  relative  strength  of  the  double-couple  force  was 
0.47  times  that  of  the  explosion. 


Amplitude  Ratios  of  the  Love  and  Rayleigh  Waves.  The  Love 
waves  generated  by  the  explosions  can  also  be  used  to  determine 
the  nature  of  the  source  mechanism.  According  to  our  formula¬ 
tion,  the  Love  waves  will  be  generated  by  the 
tectonic  component  of  the  source.  The  explosion  itself  will 
not  generate  Love  waves  under  ideal  radially  symmetric  conditions. 
In  determining  the  radiation  pattern,  we  will  again  use 
a  normalization  scheme  to  minimize  the  effects  of  propagation  paths 
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Love  v;ave  source  that  can  be  used  as  a  reference,  v;e  will 
normalize  to  Rayleigh  waves.  This  scheme  was  successfully 
applied  to  explosions  and  earthquakes  by  Toksfiz  et:  al.  (1965)  . 

The  ratio  of  the  Love  wave  amplitude  to  the  Z  compo¬ 
nent  of  the  Rayleigh  waves  generated  by  the  explosion  con 
be  written  using  equations  (1) ,  (2) ,  (J) ,  and  (4) . 


F  k/z  A,  cos  29 


t\  It:* 


(1  +  F  sin  20)kRz  Ar  (U*/Wq) 


exp ( -r  (yl 


In  writing  (7) ,  it  was  assumed  that  the  source  time  functions. 
T(w)  and  T' (u)  were  the  same  for  both  the  explosion  and  the 
double-couple  component.  Furthermore,  as  ;  the  previous  sec¬ 
tion,  it  is  assumed  that  the  tectonic  contribution  can  be 
represented  as  a  double-couple.  Thus  the  Love  waves  are 
those  generated  by  a  double-couple,  and  the  ‘Rayleigh  •.••aves 
are  the  vectorial  sum  of  those  due  to  the  explosion  and  due 
to  the  double-couple. 

In  using  (7)  we  must  compute  k^,  A^,  kR,  AR,  (U*/WQ) , 
yT  ,  Y r> •  Those  quantities  arc  functions  of  the  frequency  for  a 
given  earth  structure.  Fortunately,  k^2  Aj/kR2  AR  is  a  ~  no  tonic 
and  slowly  varying  function  of  frequency  in  the  spectral 
range  of  our  interest.  Thus  the  effect  of  the  structure  is 
minimized  lay  the  normalization  process,  and  an  imprecise 
knowledge  of  structure  does  not  limit  the  applicability  of 
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the  method.  We  took  one  average  structure  for  the  Western 
United  States  given  by  Alexander  (1963)  and  computed  the 
amplitude  response  for  the  Rayleigh  (AR)  and  the  Love  (A^) 
waves.  The  method  and  programs  of  Harkrider  (1964)  were 
used  in  these  computations.  The  results  are  shown  in 
Figures  5  and  6. 

In  computing  the  Love/Rayleigh  amplitude  ratios  from  the 
long  period  recordings  on  the  LRSM  stations  we  used  the  peak 
amplitudes.  In  cases  where  spectra  were  computed/  they  peaked 
at  about  T  =  18  seconds.  Spectral  ratios  are  shown  in  Figure 
7  for  four  stations,  together  with  the  ratio  of  peak  ampli¬ 
tudes.  On  the  average  the  agreement  is  good  enough  to  justi¬ 
fy  the  use  of  peak  amplitudes,  with  the  understanding  that 
the  results  are  valid  only  in  the  T  =  15-20  second  period 
range.  At  longer  periods  UL/URz  ratios  generally  tend  to 
increase  as  seen  for  stations  CP-CL  and  CP02  in  Figure  7. 

All  the  available-;  UL/URz  data  from  Bilby  are  shown  in 
Figure  8,  as  well  as  che  theoretical  curve  based  on  equation 
(7) .  The  reference  plane  for  double-couple  orientation 
(plane  of  0^  =  340°)  and  the  relative  strength  F  =  0.47  are 
the  same  values  that  were  determined  from  the  explosion/collapse 
Rayleigh  wave  ratios.  The  agreement  between  the  observed  data 
and  the  theoretical  curve  can  be  considered  good. 

The  consistency  of  one  single  source  model  for  both  the 
Rayleigh  and  the  Love  wave  radiation  patterns  is  encouraging. 
Although  this  is  not  a  definite  proof,  if  gives  support  to 


our  method  of  synthesizing  the  source  and  to  our  assumptions. 


11*2*  Source  Time  Function  of  Bilb^ 


The  source  time  function  of  the  explosion  car.  be  deter¬ 
mined  from  the  recordings  of  the  motion  at  distant  stations 
by  correcting  for  the  instrument  response  and  the  response 
of  the  propagation  medium.  For  surface  v/aves,  we  can  com¬ 
pute  both  amplitude  (corrected  for  attenuation)  and  phase  re¬ 
sponse  if  the  structure  is  known.  In  this  study  we  will 
use  an  average  structure  and  use  only  the  amplitude  spectra, 
since  the  accuracy  of  phase  spectra  depends  very  strongly 
on  exact  knowledge  of  the  structure. 

The  Rayleigh  waves  recorded  at  three  stations  are  first 
corrected  for  the  instrument  response  to  obtain  the  true  ground 
displacement.  These  stations  are  Kanab  (Utah) ,  Campo  (Cali¬ 
fornia)  and  Winnemucca  (Nevada) .  They  represent  excellent 
azimuthal  coverage  at  fairly  close  distances.  The  filtered 
Rayleigh  wave  pulses  from  two  of  these  are  shown  in  Figure  9 
and  their  Fourier  amplitude  spectra  in  Figure  10.  Ground 
displacement  spectra  (Figure  11)  are  obtained  by  correcting 
the  spectra  of  Figure  10  for  the  instrument  response  at  each 
station.  They  represent  the  product  of  the  source  amplitude 
spectrum  and  the  response  of  the  layered  medium  (i.e.  propa¬ 
gation  path)  to  Rayleigh  waves.  The  medium  response  includes 
the  effects  of  attenuation  and  the  source  depth. 

To  determine  the  source  function,  we  must  correct  the 
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ground  displacement  for  the  propagation  factor.  This  was 

done  using  the  impulse  response  of  the  medium  (Figure  5)  as 

given  by  equation  1.  The  effect  of  attenuation  was  removed 

using  yR  =  rf/UQ  where  U  is  the  group  velocity  and  Q  was 

assumed  to  be  IOC,  independent  of  frequency.  The  corrected 

spectra  are  shown  in  Figure  11.  This  presumably  represents  the  spectrum 

of  the  source  pressure  function. 

The  interpretation  of  amplitude  spectra  in  terms  of  a 
time  function  requires  incorporating  either 

phase  data  or  some  other  constraint.  We  - 1 

assume  that  the  pressure  pulse  has  the  form  p(t)  -  p  t) . 

This  formulation  was  discussed  in  an  earlier  study  (T«_ 
et  al . ,  1964).  The  whole  problem  now  consists  of  determining  the 
parameter  n  from  the  spectra.  A  value  of  n  =  1.5  seems  to 
agree  well  with  the  observations  as  shown  in  Figure  12.  Dis¬ 
crepancies  at  periods  longer  than  T  =  30  seconds  arc  attributed 
to  the  effects  of  lateral  heterogeneities  and  low  signal  to 
noise  ratio.  The  time  function  is  given  in  Figure  13. 
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We  must  note  here  that  the  p(t)  we  determined  represents 
the  stress  wave  form  not  at  the  source  but  at  some  distance 
from  the  point  of  detonation.  Since  wc  used  a  linear  theory 
based  on  infinitesimal  strains  to  correct  for  propagation  and 
attenuation  effects,  our  corrections  arc  valid  only  to  the 
boundary  of  the  region  where  these  conditions  are  met.  This 
may  be  at  a  distance  of  Several  kilometers  from  the  source- 
point.  We  must  clarify  one  other  aspect  of  Figure  9:  the 
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pressure  pulse  was  based  on  data  in  the  period  range  of  10  to 


35  seconds,  so  we  coxild  not  see  any  of  the  fine  features  of 


the  pulse  that  would  be  observable  primarily  in  the  high 


freouenc''  components.  Neither  could  ve  see  a  small  residual 


cermanent  displacement  or  strain.  Even  with  these  limita¬ 


tions  the  share  of  the  time  function  shown  in  Figure  9  is 


similar  to  (although  somewhat  broader  than)  those  of  close-in 


measurements  (Wistor,  et  al.,  1963;  Perrett,  1963).  The  source 


time  function  varies  since  it  is  dependent  on  the  yield  of  the 


explosion,  the  medium,  and  shot  depth.  In  gene-  larger  ex¬ 


plosions  have  broader  pulses  (smaller  n  values)  r  investi¬ 


gators  have  also  proposed  a  step-function  for  the  source  model 


on  the  basis  of  surface  wave  data  (Tsai  and  Aki,  1970).  However, 


the  available  close-in  acceleration  measurements  (Perrett,  1968) 


and  near-field  strain  data  (Smith,  et  al. ,  1969)  indicate  larger 


stresses  at  the  beginning  of  the  pulse  than  toward  the  tail. 


These  emphasize  the  need  for  spectral  data  near  the  source  over  a 


broad  period  range  for  future  detailed  source  studies, 


It  would  be  of  .interest  to  examine  the  source  function 


of  the  Love  waves  generated  by  Bilby.  Choosing  four  stations 


where  Love  ’waves  are  well  separated  from  Rayleigh  wave  inter¬ 


ference,  wo  followed  the  same  procedure  as  for  Rayleigh  waves 


of  correcting  the  spectra  for  the  instrument  response,  and 


the  response  of  the  layered  medium  for  an  orthogonal  doublc- 


couole  source.  ihe  resultant  spectra  of  the  source  time 


function  are  shewn  in  Figure  1.4,  and  their  shapes  are  con¬ 


sistent  for  all  four  stations. 


From  the  comparison  of  Figures  12  and  14  it  is  obvious 
that  the  source  time  function  of  the  explosion-generated 
Rayleigh  v/aves  and  that  of  the  Love  waves  is  quite  different. 

The  Love  v/aves  seem  to  be  richer  in  low  frequency  components. 

From  the  available  amplitude  data  in  a  relatively  narrow  fre¬ 
quency  band,  v/e  cannot  determine  the  time  function  of  Love 
waves.  It  appears  to  be  between  a  step  function  and  a  ramp, 
if  v/e  assume  a  point  source  near  the  surface,  a  value  of  Q  =  TOO, 
and  if  v/e  ignore  the  effects  of  source  volume.  With  these  limi¬ 
tations,  a  step  function  w*  th  a  rise  time  on  the  order  of  a  few 
seconds  may  be  an  acceptable  approximation  for  Love  v/aves. 

The  source  spectra  of  the  Rayleigh  and  Love  waves  generated 
by  explosions  are  important  because  they  contain  information 
about  the  mechanism  of  the  Love  wave  generation.  The  differ¬ 
ences  demonstrated  above  are  significant.  They  indicate  that 
an  explosion  may  trigger  the  Love  wave  radiation,  but  it 
probably  does  not  control  its  time  history.  The  Love  wave 
source  spectra  are  closer  to  those  of  earthquakes  tnan  to  the 
spev'tra  of  explosion  generated  Rayleigh  waves. 

I I . 3 .  Comparison  with  other  Explosions 

In  determining  the  source  properties  of  the  Bilby  explo¬ 
sion,  we  followed  fc»;e  same  procedures  that  were  used  in  our 
earlier  studies  of  the  Hardhat,  Haymaker,  Sedan  and  Shoal 
explosions,  as  well  as  in  some  later  events  (Toksdz,  et  ai., 
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1965;  Kehrer,  1969;  Toksdx  et_  al_.  ,  197|).  Here  we  will  compare 
the  Bilby  results  with  those  of  other  explosions,  which  are 
tabulated  in  Table  1. 

The  most  significant  result  is  that  Bilby,  like  the 
Haymaker  and  Shoal  explosions ,  generated  Love  waves .  The 
source  mechanism  in  all  cases  can  be  explained  in  terms  of  an 
ideal  radial  explosive  source  superimposed  on  a  tectonic  source 
of  double-couple  form.  The  orientation  and  relative  strength 
of  the  double-couple  seem  to  be  controlled  by  the  properties  of 
the  medium  and  the  orientation  of  the  tectonic  axes  in  the  source 
region.  Bilby  fin  tuff)  and  Haymaker  (in  alluvium)  were  located 
about  5  km  apart.  The  radiation  patterns  of  Rayleigh  waves  are 
almost  identical.  In  both  cases  the  principal  plane  of  the 
double-couple  is  oriented  in  the  direction  0  =  340°.  Other 
explosions  in  the  same  general  area  give  similar  orientation. 

In  the  shoal  explosion,  which  was  fired  in  a  completely  differ¬ 
ent  area,  the  orientation  of  the  double-couple  was  in  very  good 
agreement  with  those  of  earthquakes  in  the  area  (Toksfiz  et  a_L .  , 
1965) .  These  facts  suggest  that  multipolar  contributions  to 
the  radiation  patterns  arc  controlled  by  the  general  tectonic 
features  of  the  region. 

The  relative  strength  of  the  tectonic  (double-couple) 
contribution  to  the  radiation  pattern  seems  to  he  controlled 
by  the  properties  of  the  medium  in  which  the  explosion  is  de¬ 
tonated.  This  is  demonstrated  by  the  observation  that  for  ex¬ 
plosions  in  salt  domes  (Gnome  and  Salmon)  and  loose  alluvium 
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(Sedan) ,  the  source  functions  did  not  have  multipolar  compo¬ 
nents  (i.e.,  F  -  0).  For  Haymaker,  v/hich  was  buried  deep  in 
alluvium,  a  value  F  =  0.33  was  determined.  For  Bilfoy  (in 
tuff)  F  was  0.47,  and  for  Shoal,  which  was  fired  in  granite, 

F  was  equal  to  0.9.  Hardhat  was  another  explosion  in  granite, 
and  it  had  an  F  value  greater  than  1.0.  These  and  t.en  other 
explosions  recently  studied  (Toksdz,  o_t  al .  ,  lS7,rj  indicate 
an  increase  of  F  with  increasing  rook  strength  and,  hence,  strain 
energy  capacity  of  the  medium. 

From  these  examples  we  cannot  determine  conclusively 
v/hether  the  multipolar  component  of  the  seismic  energy  radia¬ 
tion  is  due-  to  release  of  some  of  the  strain  energy  accumu¬ 
lated  in  the  medium,  due  to  fracturing  or  due  to  earthquake 
triggering.  Our  laboratory  model  experiments, 
described  in  the  next  section, will  help  to  clarify  some  of 
these  points. 
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I  III.  LABORATORY  EXPERIMENTS 
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Laboratory  experiments  can,  in  principle,  provide  a  better 
physical  understanding  of  explosive  seismic  sources  than  field 
observations  because  the  explosive  source  strength,  pre-stress 
field  and  elastic  properties  of  the  Medium  can  be  varied  in  a 
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known  way.  Two  separate  laboratory  experiments  were  carried 
out  with  explosions  in  pre-s tressed  media.  The  first  experi¬ 
ment  utilized  a  rod  stressed  in  torsion,  an  explosive  source, 
and  strain  gauges.  The  second  set  of  experiments  was  per¬ 
formed  with  point  sources  in  pre-stressod  two-dimensional 
plate  model.  Photoelastic  techniques  and  high  speed  photo¬ 
graphy  were  utilized  to  determine  the  shot  effects  in  the 
source  region,  radiation  of  stress  waves,  and  the  crack  for-  ■ 
mation  as  a  function  of  medium  properties  and  pre-stress 
levels. 

III.l.  Stressed  Rod  Experiment 

A  simple  one-dimensional  laboratory  experiment,  giving 
insight  into  the  basic  physical  mechanisms  involved  in  the 
production  of  stress  waves  from  explosive  sources  in  pre¬ 
stressed  media,  is  shown  in  Figure  15.  Plexiglas  rods  1.27 
cm  in  diameter,  and  of  varying  lengths  up  to  2  m,  were  stressed 
in  torsion.  A  few  cm  from  one  end,  a  bridge  wire  was  inserted 
in  a  close  fitting  diametric  hole  and  exploded  with  a  100- 
joule  source.  This  was  adequate  energy  to  sever  the  end  of 
the  rod  leaving  the  remainder  of  the  rod  free  from  the  torsional 
restraints  holding  it.  The  resulting  radiated  waves  were 
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recorded  by  at''ni.n  gauges  placed  on  the  rod  cireumf cronce 
at  various  distances  from  the  shot  but  oriented  at  45°  to  the 
rod  axis.  Typically  four  gauges  were  used,  one  recording 
directly  onto  an  oscilloscope  and  all  four  were  recording  on 
magnetic  tape.  The  gauges  were  also  read  statically  before 
firing  the  bridge  wire  source  to  provide  a  measure  of  the 
prestrain. 

Typical  oscillogram  tracings  from  four  shots  at  dif¬ 
ferent  pre-stress  levels  are  shown  in  Figure  16.  Two  out¬ 
standing  pulses  on  the  oscillogram  are  the  one-dimensional 
equivalents  of  P  and  S  waves.  Velocities  measured,  using 
gauges  at  several  distances  from  the  source,  identify  these 
as  the  rod  wave  (P)  and  torsion  wave  (S) ,  respectively.  The 
distinctive  wave  shapes  for  P  and  S  pulses  in  Figure  16  in¬ 
dicate  that  P-wave  shape  is  controlled  by  the  time  function 
of  the  explosive  lourca .  The  S-wave  is  characterized  by  a 
step-function  input,  and  corresponding  to  the  permanent  re¬ 
laxation  of  the  torsional  strain  in  the  medium. 

A  number  of  oscillograms  were  obtained  at  various  pre¬ 
strain  levels  but  with  the  same  explosive  source.  The  ampli¬ 
tudes  of  the  P  and  S  waves  are  plotted  as  a  function  of  the 
pre-strain  level  in  Figure  17.  It  is  quite  apparent  that  the 


amplitude  of  the  S-wave  duo  to  strain  release  increases 
linearly  with  an  increase  in  the  prestressing  field  while 
the  P-wave  amplitude  remains  constant.  From  the  zero  inter¬ 
cept  times  for  both  P  and  S  waves  in  time-distance  plots,  it 
can  be  concluded  that  at  the  precision  level  of  this  experiment. 
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there  was  no  evidence  of  a  different  origin  time  for  P  than 
for  5.  Furthermore ,  the  observed  time  function  of  radiated 
S-waves  {Fig.  16)  can  bo  characterised  as  a  Heaviside  step 
function,  as  seen  through  response  of  strain  gauges.  This 
is  consistent  with  the  source  time  functions  observations 
for  Lhe  Rayleigh  waves  and  the  Love  waves  from  the  Bilby 
exp]  osier.,  as  determined  in  tho  previous  section. 


III. 2.  Two-dimensional  Photoelastic  Experiments 


The  one-dimensional  experiment  cannot  be  used  to  model 
the  spatial  variation  of  radiation  pattern  from  explosive 
sources  in  pre-s tressed  media.  However,  a  close  simulation 
of  the  underground  explosions  can  be  achieved  by  using  two- 
dimensional  modeling  techniques. 

The  experimental  procedure  used  in  this  study  involves 
glass  and  plexiglas  plates  stressed  in  tension  or  shear.  The 
explosive  source  is  simulated  by  a  small  radial  shot  and  the 
source  and  radiation  fields  are  observed  by  dynamic  photo¬ 
elastic  techniques.  This  experiment  differs  from  previous 
studies  by  Kim  and  Kissiinger  (1967)  utilizing  aluminum  and 
plexiglass 
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Mia*'*.::  iivl  .  :  >  .in  g-iuge  techniques.  The  *  ii  ui f ferences  are 
that  in  the  current  study:  1.  A  pre-stressi ng  frame 
capable  of  synthesizing  a  variety  of  static  stress  fields 
including  tension,  compression,  and  shear  '/as  used.  2.  The  details 
of  the  explosion,  dynamic  cracking,  and  wave  generation  in 
the  near  source  region  were  continuously  observed  through 
the  use  of  dynamic  photoeiastic  techniques.  3.  Glass  was 
included  among  the  model  material  used.  The  shock  loading 
properties  of  glass  are  closer  to  those  of  rock-forming 
silicates,  and  its  more  brittle  character  simulates  near 
surface  rocks  somewhat  better  than  the  more  ductile,  but 
.experimentally  more  convenient,  plastics  cr  metals.  In  the 
experiment,  glass  turned  out  to  be  more  important  than 
Plexiglas  for  studying  cracking  phenomena  because  the  c rack 
velocity  in  glass  was  higher  than  the  expansion  velocity  of 
explosion  gases,  whereas  in  the  Plexiglas  it  was  lower.  There¬ 
fore,  only  in  the  case  of  glass  did  we  obtain  adequate  de¬ 
tailed  photographic  evidence  of  tha  dynamic  cracking  phenomena. 

Figure  i8  shows  in  outline  the  apparatus  used  in  our  two- 
dimensional  work.  Referring  to  this  figure,  explosive  excita¬ 
tion  of  the  model  was  provided  by  a  closed  loop  of  mild  detonating 
fuse  (MDF) .  This  induced  a  radial  shock  source  in  the  plate 
as  the  shock  wave  from  the  detonation  in  the  MDF  enetered  the  plate 
model  simultaneously  from  both  sides.  The  model  was  surrounded  by 
polarizing  optics  so  that  the  e.. plosion  process,  cracking, 
wave  generation  and  propagation  could  be  made  visible  to  r.he 


high  spaed  '.'.Ming  camera  by  means  of  the  photoelastic  effect 
Suitably  conditioned  intense  paraxial  illumination  was  pro¬ 
vided  by  the  Argon  bomb,  a  front  surface  mirror,  and  the 
Fresnel  lens. 

A  much  more  detailed  discussion  of  the  experimental 
technique,  the  methods  of  using  dynamic  photoelastic  obser¬ 
vations  (isochromatics  and  isoclinics)  to  identify  seismi- 
cally  interesting  features,  is  given  in  our  earlier  paper 
(Thomson  sit  al_.  ,  1969).  and  will  not  be  further  elaborated 
on  here. 

In  addition  to  photoelastic  observations,  strain  gauges 
were  used  to  record  the  dynamic  field  at  points  too  far 
distant  from  the  source  for  the  optical  techniques  to  be 
completely  adequate.  These  gauges  were  set  out  so  as  to 
instrument  one  quadrant  surrounding  the  source. 

In  the  following  section  we  will  investigate  in  detail 
both  the  propagation  of  the  stress  waves  from  the  shot  and 
the  characteristics  of  the  cracking  phenomena. 

III. 2.1.  Stress-Wave  Propagation  in  a  Stressed  Plate 

To  understand  the  fringe  patterns  shown  in  Figures  19 
to  24,  it  is  necessary  to  introduce  some  basic  definitions 
and  results  of  some  calculations. 

Upon  passage  through  a  stressed  plate  model,  monochro¬ 
matic  circularly  polarized  light  becomes  polarized  in  two 
orthogonal  directions. 


These  are  parallel  to  the  principal 
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in  l  ho  plane  of."  the  p'laln.  ;‘.i nou  the  Light 
polarized  in  each  direction  experiences  a  different  velocity 
in  traversing  the  model,  an  interference  effect  due  to  the 
resulting  phase  difference  is  observed  v/hen  the  two  compo¬ 
nents  of  the  incident  beams  are  combined  upon  passage  through 
the  analyzer.  The  fringes  so  produced  are  called  isochroma¬ 
tics  and  are  proportional  to  the  difference  in  the  principal 
stresses  (Oj  -  o 2)  and  hence  the  maximum  shearing  stress. 

The  fringe  order  is  thus  given  by: 


n  - 


h  (o j  -  a2)/f 


(8) 


n  - 


2h 

f 


jf  | 


where  t  is  the  maximum  shear  stress,  h  is  the  plate  thickness 
and  c  is  the  constant  of  proportionality  which  is  wavelength 
dependent.  Therefor^  if  white  light  is  used  to  illuminate 
the  model  the  isochromatic  effect  is  a  multicolored  inter¬ 
ference  pattern.  Knowing  the  stress  field  {or,  oQ,  arQ)  one 
can  calculate  the  isochromatics  from  equation  (8)  and  to 
within  a  numerical  constant  factor  (for  a  particular  material 
and  thickness) .  This  factor  must  bo  determined  experimentally. 

In  plane  polarized  light  not  only  arc  the  isochromatics 
observed  but  an  additional  set  of  extinctions,  called  the 
isoclinics,  are  seen  wherever  the  principal  stress  directions 
are  parallel  or  perpendicular  to  the  pass  directions  of  the 
polarizer  and  analyzer.  In  this  experiment,  polarizer  and 
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n'od  at  45°  anJ  135°. 


Thu?;  in  the  unstressed 


v.'.iv-  o •' ■  • 

plate  case,  shot  generated  P  waves  will  have  isoclinics  in 
directions  S  -  45,  135,  225,  315°.  The  S-waves  from  the  same 
source  will  have  isoclinics  at  0  =  0,  90,  180  and  270°.  The 
isoclinics  are  used  to  identify  the  wave  types  (P  or  S)  in 
these  experiments. 

In  the  case  of  the  stressed  plate  the  isoclinics  will 
shift  due  to  the  static  stress  field.  In  this  case  the  loci 
of  isoclinics  can  be  calculated  from  the  vectorial  summation 
of  principal  static  and  dynamic  stresses. 

Examples  of  isochromatics  and  isoclinics  from  a  shot  in 
.an  unstressed  plexiglas  plate  are  shown  in  Figure  19.  Note 
that  these  patterns  are  very  clear  for  P-waves,  but  less 
clear  for  S-waves  which  were  much  smaller  in  amplitude.  In 
fact,  S-waves  can  be  identified  in  this  case  on  the  basis  of 
the  isoclinics.  The  theoretical  isochromatic  and  isoclinic 
pattern  for  a  radial  transient  source  in  a  plate  xs  given  in 
Figure  20. 

In  performing  the  calculations  on  which  Figure  20  (and 
also  Figure  23}  wore  based,  it  was  assumed  that  the  dynamic 
stress  components  displayed  geometric  attenuation  inversely 
proportional  to  the  square  root  of  the  distance  from  the 
source  and  that  the  wave  form  propagated  without  shape  change. 
These  assumptions  are  only  approximations  and  in  order  to 
better  understand  the  effect  on  the  i  sochromatics  v/e  have 
made  a  few  precise  theoretical  calculations  for  the  cylindri¬ 


cal 


source  in  a  plate.  The  analytic  formulations  of 


chis  problem  have  been  given  by  Kromm  (1948).  The  numerical 
solution  of  the  integral  equation  has  been  discussed  by 
Barker/  Toksfiz,  and  Ward  (1959).’  The  results  for  a  particu¬ 
lar  transient  cavity  pressure  are  shown  in  Figure  21,  where 
the  radial  stress,  displacement,  and  maximum  shear  stress 
have  been  computed  as  a  function  of  normalized  distance 
r  =  distance/hole  radius  from  the  source  for  glass  and 
plexiglas.  Looking  at  the  maximum  shear  stress  in  Fig.  21, 
it  is  clear  that  fringes  would  separate  better  in  the  case 
of  plexiglas  than  in  the  case  of  glass. 

The  principal  features  of  Figure  19  were  identified  by 
comparing  them  to  the  theoretical  isochromatics  and  isoclinics 
of  Fig.  20  and  by  the  use  of  time-distance  plots  to  obtain 
velocity.  Fractional  fringe  orders  only  (n  <  1)  are  visible 
in  this  framing  camera  sequence.  The  compressive  and  dila- 
tational  portions  of  P,  the  zero  order  isochromatic  separating 
them  and  a  weak  S  waves  are  all  visible  in  this  photograph. 

The  spreading  effect  on  the  isoclinics  caused  by  superposition 
of  P  and  S  waves  indicated  theoretically  in  Fig.  20  can  be 
seen  in  fact  by  close  inspection  of  Fig.  19. 

Using  black  and  white  photography,  the  actual  stress 
level  is  determined  from  the  observed  fringe  order  n  and 
equation  (8) .  Color  photography  was  used  to  compute  stress 
levels  from  the  retardation  of  one  wavelength  with  respect, 
to  another  wavelength  and  the  resulting  chromatic  effects. 
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The  main  emphasis  of  this  s/ork  is  .to  determine  the  effect 
cf  pre-stress  on  the  generation  and  radiation  of  seismic  waves. 
For  this  reason,  photoelastic  experiments  with  explosive 
sources  were  repeated  with  prestressed  plates.  The  plexi- 
glas  and  glass  plates  were  stressed  in  pure  tension  or  pure 
shear  using  the  stress  frame.  The  shear  was  generated  by 
applying  equal  amounts  of  tension  and  compression  in  the  y 
and  x  directions,  respectively.  The  maximum  amount  of  stress 
was  limited  to  about  100  bars  per  axis. 

The  static  stress  field  is  somewhat  complicated  by  the 
presence  of  a  small  hole  {drilled  to  place  the  explosive 
source)  in  the  center  of  the  plate.  The  radial,  cr  ,  tangen¬ 
tial,  a0  ,  and  shear,  Tr0  ,  components  of  the  stress  around 

the  hole  of  radius,  a,  can  be  computed  from  Muskhelishvili's 
* 

solution  (Savin,  1966). 

In  our  experiments  the  hole  radius  is  a  =  0.185  cm.  Thus 
in  most  of  the  plate  where  fringe  patterns  are  observed,  p  <<  1. 
Neglecting  the  higher  order  terms  of  p  ,  the  static  stress 
field  can  be  expressed  as  (Thomson,  et  al. ,  1969) 


a  -  —  (1  +  cos  20) 
r  2 

aQ  =  j  (1  -  cos  29) 


T 


r0 


a 

—  sm 
2 


29 


(9) 
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for  tension,  and  us 

at  ~  v  cos  28 

a  -  -a  cos  2<3  (10) 

8 

t j.0  =  -o  sin  20 

for  simple  shear  stress,  where  8  is  measured  counter-clockwise  * 
from  tensile  axis. 

When  an  explosive  source  is  fired  in  the  pre-stressed 
plate,  the  dynamic  stress  components  are  superimposed  over 
the  static  field  given  above,  producing  a  complicated  iso- 
chromatic  and  isoclinic  pattern  (Thomson,  et.  al . ,  1969). 

These  complications  make  it  difficult  to  extract  the  dynamic 
radiation  patterns  from  the  framing  camera  photographs,  and 
therefore  strain  gauge  observations  were  used  for  this  purpose. 
The  radiation  pattern  of  seismic  waves  from  an  explosive  source 
in  a  glass  plate  under  tensile  stress  is  shown  in  Fig.  22. 

In  this  figure,  ir.  addition  to  a  complicated  S-wave  radiation 
pattern,  we  also  see  the  development  and  growth  of  cracking 
and  radiation  from  crack  tips.  These  experiments  were  repeated 
both  with  plexiglas  and  glass  plates  under  shear  pre-stress 
conditions.  The  radiation  patterns  were  similar  to  the  above 
cases.  In  the  case  of  pre-stressed  plates,  the  shear  waves 
were  more  abundant  than  in  the  unstressed  plate  experiments. 

The  maximum  stress  level  was  determined  by  the  static 
breaking  strength  of  the  plate,  which  was  lowered  considerably 
because  of  the  presence  of  the  shot  hole.  Furthermore,  in  the 
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case  of  glass  p.latos  static  flexural  failure  occurred  before 
a  tensile  or  shear  failure.  In  the  dynamic  experiments ,  the 
maximum  static  stress  level  achieved  was  167  bars. 

III. 2. 2.  Strain  Guage  Measurements  of  P-S  Wave  Radiation  Patterns. 

A  set  of  four  BLH  foil  strain  gauges  were  used  t6  determine 
the  azimuthal  variations  of  P  and  S  wave  radiation  patterns  under 
different  pre-stress  conditions.  Gauge  outputs  were  recorded  on 
oscilloscopes  writing  in  the  single-sweep  mode.  Each  gauge  had  an 
active  area  of  8x10  mm,  and  this  limited  its  response  to  periods 
greater  than  4  microseconds.  In  most  of  the  experiments  strain 
gauges  were  oriented  at  45°  to  the  radial  direction  and  placed 
at  5  azimuths,  at  a  =  0,  30,  60,  and  90°,  at  a  distance  of 
16  cm  from  the  source.  Since  the  presence  of  symmetry  was 
shown  from  framing  camera  pictures,  instrumenting  only  one 
quadrant  with  strain  gauges  was  sufficient  for  determining 
the  radiation  pattern. 

The  strain  gauge  experiments  we re  carried  out  with  shots 
in  both  unstressed  and  stressed  plates.  In  the  case  of  un¬ 
stressed  plates,  no  significant  S-waves  wore  detected  above 
the  noise  level.  In  the  case  of  pre-stressed  plexiglas  and 
glass  plates,  however,  well-defined  S-waves  were  recorded. 

The  radiation  patterns  of  P  and  S-waves  for  plexiglas  and 
glass  plates,  stressed  under  tension  to  163  bars  and  114  bars, 
respectively,  are  shown  in  Figure  23. 

Utilizing  the  local  plane  wave  formulation,  strain  energy 
can  be  estimated  from  the  observed  compress ional  and  shear 
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wave  strains  radiated. 


3 

W  =  2  p.e.  (11) 

i=l  •L  * 


N 

where  W  =  energy  density  and  p.  and  e.  are  principal 

stresses  and  strains  respectively.  For  compressional  waves 

the  strain  measured  with  a  gauge  oriented  at  45°  to  the 

radial  direction  (e)  is  related  as  =  2e  and  tor  shear 

waves  e  =  =  ~e2*  Thus  the  compressional  and  shear  wave 

potential  energy  density  fluxes ,  W  and  W  ,  can  be  expressed 

p  s 

as 


V  —■ — (AlkL.  e2  =  H  e2  V, 
p  X+2p 


(12) 


W„  =*  2  p  e2  V 

s  s 


where  V  and  V  are  plate  and  shear  velocities,  respectively, 
p  s 

For  plexiglas  and  glass  used  in  this  experiment,  M  =  0.95  x  10 
11  2 

and  13.2  x  10  dynes/cm  ,  respectively.  The  corresponding 

rigidities  for  plexiglas  and  glass  are  y  =  0.14  3  x  10^  and 
11  2 

2.5  x  10  dynes/cm  .  The  shear- to -compressional  potential 
energy  density  ratios  are  shown  in  Figure  24  for  peak  values 
’  of  e . 

From  the  radiation  patterns  of  both  the  peak  energy 
density  ratios  and  the  strain  amplitudes  themselves  (Figures 
23  and  24),  it  is  clear  that  the  S  wave  radiation  can  be 
explained  by  either  a  dipole  or  double  couple  type  source 
having  approximately  a  sin  20  azimuthal  dependence.  The  P 


to 
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wave  radiation  is  less  clear  since  the  contribution  from  the 
explosion  as  well  as  radiation  due  to  pre-stress  are  super¬ 
imposed.  It  is  apparent  that  the  radiation  pattern  for  P 

would  be  approximated  by  a  function  of  8  of  the  form 
2 

(A  +  B  cos  0)  where  A  and  B  are  constants,  suggesting  again 
the  dipole  source  with  the  addition  of  an  isotropic  component. 

Inspecting  the  radiation  patterns  and  energy  ratios 
displayed  in  Figures  23  and  24,  we  are  faced  with  two  questions! 
(1)  where  does  the  S  wave  energy  come  from,  and  (2)  what  is  the 
reason  for  the  skewness  in  the  P  wave  radiation  patterns? 

Let  us  first  consider  the  energy  question.  The  total  energy 
of  an  elastic  wave  is  equal  to  twice  the  potential  strain  energy 
integrated  over  a  cycle.  Thus  we  can  compute  the  P  and  S 
wave  energies  by  integrating  the  energy  fluxes  given  in  (12) 
for  the  duration  of  the  pulse  first,  and  then  integrating  over 
the  radiation  pattern. 

At  the  level  of  approximation  being  used  here  we  can 
evaluate  these  integrals  by  fitting  each  of  the  compressional 
and  dilatational  half-cycles  of  strain  wave  forms  with  a  half 
sine  wave.  Integrating  the  flux  azimuthally,  allowing  for 
the  radiation  pattern,  we  obtain  approximate  expressions  for 
the  total  energy  E  and  E  in  P  and  S  waves  for  each  half- 

r  w 

cycle  of  duration  t.g : 

Ep  =  »hpV>  Aj5  t0  r 

(13) 

Es  -  »hpv^  Bjs  t0  r 
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where  h  is  the  moc'ul  thickness,  r  is  the  distance  of  the 
strain  gauqe  from  the  source  and  A^g  end  B^  are  peak  strain 
gauge  readings  for  the  P  and  S  waves  respectively  in  the 
record  interval  under  study,  and  the  gauges  are  assumed 
located  at  the  peak  of  the  radiation  pattern  and  oriented  at 
tt/4  to  the  radius  vector  from  the  source.  For  gauges  not 
located  at  the  peak  of  the  radiation  pattern  an  additional 
correction  was  applied  in  obtaining  total  energy  values  as 
discussed  subsequently  in  connection  with  the  radiation 
patterns.  For  three  shots  shown  in  Figures  23  and  24  (shots 
No,  13028,  13029,  13030),  the  corresponding  P  wave  energies 
(averaged  from  the  observations  at  each  gauge)  are:  313,  231 
and  143  x  10  ergs.  The  S-wave  energies  are  49.3,  3.72  and 
0.73  x  10^  ergs.  The  corresponding  tensile  prestress  figures 
in  bars  are  163,  1S3,  and  114.  Note  that  the  first  two  shots 
are  in  plexiglas  and  the  third  in  pyrex  glass. 

The  S-wave  energies  cannot  come  from  the  conversion  of 
the  explosive  energy  into  S  waves  since  very  little  if  any 
S  waves  are  generated  by  the  explosions  in  unstressed  plates. 
Thus,  we  must  look  .into  the  energy  release  from  the  stressed 
plate.  The  explosion  crushes  a  small  circular  region  of  the 
plate  around  the  shot  point.  This  normally  does  not  exceed 
a  radius  of  1  cm.  Beyond  this  region  some  radial  cracks  extend 
to  greater  distances.  If  all  the  strain  energy  in  a  zona  of 
radius  a^  is  radiated  as  seismic  energy,  the  total  amount 
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where  E  is  Youngs  modulus  and  was  determined  from  observed 
plate  wave  velocities  and  h  is  the  plate  thickness.  Equation 
(14)  gives  a  good  approximation  to  the  upper  limit  of  the 
available  seismic  energy  if  a^  is  taken  to  be  the  radius  of 
the  "quasi-elastic"  cracked  zone  (Archambeau  and  Sammis, 

1970) .  The  problem  that  arises  in  estimating  the  energy 
from  (14)  is  that  of  the  definition  of  the  cracked  zone 
radius.  In  our  experiments,  the  approximate  values  for  the 
radius  of  the  totally  crushed  and  extensively  cracked  zone 
were  about  0.4  cm  and  1  cm,  respectively.  The  values  of  EST 
for  the  strain  levels  of  shots  13028,  13029  and  13030  are 
tabulated  in  Table  2  for  the  larger  value  of  a.  Comparing 
these  figures  with  the  previously  given  wave  energies,  it 
will  be  apparent  that  the  available  strain  energy  is  smaller 
than  the  total  energy  radiated  as  S  waves  for  a  =  1  cm. 


We  should  further  consider  that  not  all  the  available 
strain  energy  is  radiated  as  seismic  energy.  In  fact,  only 
a  small  fraction  is  likely  to  be  radiated  in  the  spectral 
range  of  our  measurements,  and  this  too  must  be  divided  into 
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P  and  S  waves.  Thus  it  is  clear  that  a  source  region  much 
larger  than  the  crushed  zone  must  be  considered.  This  leads 
us  to  consider  relaxation  of  the  medium  along  the  cracks  that 
radiate  to  greater  distances  from  the  source  than  those  con¬ 
sidered  in  Table  2.  The  cracking  phenomenon  will  be  dis¬ 
cussed  in  the  next  section. 

The  second  problem  we  must  consider  is  the  skewness  of 
the  P  wave  radiation  pattern.  This  phenomenon,  apparent  on 
our  radiation  patterns,  can  also  be  seen  on  the  data  of  Kim 
and  Kisslinger  (1967) .  It  cannot  be  attributed  to  asymmetry 
of  the  shot  since  it  occurs  only  in  the  pre-stressed  experi¬ 
ments.  Since  there  is  a  contribution  to  P  waves  from  the 
radiation  of  the  strain  energy  due  to  the  relaxation  of  the 
plate,  then  the  mechanism  of  this  relaxation  must  be  taken 
into  account.  From  the  above  discussion  as  well  as  from  the 
description  of  the  cracking  in  the  next  section,  the  growth 
of  the  cracks  does  affect  the  radiation  pattern.  Furthermore, 
where  cracks  extend  close  to  the  observation  position,  the 
growth  can  be  characterized  as  a  unilateral  tensile  fault. 

This  presumably  is  due  to  the  scattering  effect  of  the  shattered 
region  behind  the  crack  tips,  on  radiation  not  projected 


broadly  in  the  direction  of  crack  propagation.  '.■.’e  have  computed 
the  radiation  pattern  of  a  unilateral  tensile  fault  propaga¬ 
ting  perpendicular  to  the  tension  axis,  utilizing  Ang  and 
’Williams  *  (1959)  formulation  as  given  by  Savage  and  Mans  inha 


(1963)  .  The  radiation  patterns  for  P  and  S  -waves  are 
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proport ional  to: 


and 


c  {1-2V^  eos^  O/Vp)  /  (Vp  -  c  ccs  0) 


c  sin  20  /  (V-  -  c  cos  0) 


(15) 


respectively ,  where  c  is  the  crack  velocity.  In  our  experi¬ 
ments,  the  observed  radiation  patterns  of  Figure  23  and  24 
have  an  additional  isotropic  component  due  to  the  explosive 
source  affecting  the  P-wave  radiation  patterns.  A  set  of 
theoretical  radiation  patterns  were  calculated  using  (15) 
with  the  addition  of  various  isotropic  components.  These  are 
shown  in  Figure  25.  The  data  of  Figure  25  apply  to  plexi- 
glas. 

Two  aspects  of  Figure  25  are  noteworthy.  (a)  Increasing 
crack  propagation  velocity  produces  a  shift  of  the  radiation 
pattern  maximum  in  the  direction  of  crack  propagation  and 
(b)  increasing  the  isotropic  component,  while  broadening 
the  radiation  pattern,  dees  not  change  the  position  of  the 
maximum.  The  radiation  pattern  maxima  for  P  and  S  and  for 
glass  and  plexiglas  were  calculated  as  a  function  of  crack 
velocity.  The  results  are  shown  in  Figure  26. 

Inspecting  the  shift  in  the  radiation  pattern  of  P  waves 
from  the  observations  of  Figure  2  3  and  comparing  this  to  the 
theoretical  shift  of  Figure  26,  it  is  apparent  that  only  a 
crude  estimate  of  the  crack  velocity  can  be  obtained,  because 
of  the  limited  number  of  data  eh  ints  available.  A  crack 
velocity  of  six  tenths  the  shear  wave  speed  is  consistent  with 
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observations  foi:  both  materials.  The  crack  speed  was,  however. 


measured  accurately  from  the  optical  observations,  and  the 
subject  of  crack  propagation  will  be  dealt  with  separately  in 
the  next  section.  In  reducing  the  strain  gauge  observations 
to  wave  energy,  as  listed  in  Table  2,  gauge  observations 
wore  corrected  for  the  radiation  pattern  effects  using  a 
crack  velocity  0.6  Vg, 


IT. 2. 3.  Near  Source  Phenomenon  and  Crack  Propagation 

The  stresses  in  the-  immediate  vicinity  of  the  explosion 
are  very  high  and  may  reach  hundreds  of  kiiobars.  This  easily 
exceeds  the  compressive  strength  of  the  material,  and  a  crushed 
zone  which  extends  from  a  few  millimeters  to  a  centimeter  or 
more,  depending  on  the  shot  size  and  the  medium,  is  developed. 

As  the  stress  wave  propagates  outward  from  the  shot  point, 
the  pressure  falls  gradually  until  the  elastic  limit  of  the 
material  is  reached.  Between  the  crushed  and  elastic  zones, 
there  is  a  region  of  radial  cracking  produced  by  the  tensile 
failure  during  the  unloading  (dilatationai)  phase  of  the  stress 
waves.  Because  of  the  importance  of  this  region  for  the 
radiation  patterns  of  seismic  waves,  great  emphasis  was 
placed  on  understanding  its  properties  in  our  dynamic  photo- 
elastic  experiments. 

* 

The  immediate  source  region  of  an  explosion  xn  g  "  thick 
glass  plate  is  shown  in  Figure  27.  The  evaporated  zone  (totally 
dark)  and  the  crushed  zone  are  indicated  very  clearly.  The 
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ini  ‘.vial  radial  cracks  are  visible  along  the  inner  edges  of  the 
isochromatics . 

In  unstressed  plate  experiments  the  radial  cracks  are 
distributed  isotropically  and  grow  in  radial  directions.  In 
the  case  of  pre-stressed  plates#  initially  the  cracks  do 
start  radially.  However#  after  some  time  the  cracks  either 
stop  growing  or  turn  in  the  direction  dictated  by  the  stress 
axes.  In  the  case  of  plates  under  tension,  the  favored  direc¬ 
tion  of  growth  is  perpendicular  to  the  tension  axis. 

The  propagation  of ' explosion-induced  cracks  in  glass 
plate  under  tension  is  shown  in  Figure  28.  The  preferred 
growth  (in  direction  perpendicular  to  stress  axis)  and  the 
radiation  of  seismic  energy  from  the  cracks  (light  zones)  are 
indicated  clearly.  The  stress  concentration  at  crack  tips 
is  visible.  The  dark  areas  are  the  stress -free  (relaxed)  re¬ 
gions  . 

In  addition  to  the  above  observations  we  have  been  able 
to  determine  the  velocity  of  the  crack  propagation  by  follow¬ 
ing  the  progress  of  crack  tips  in  successive  frames.  By  re¬ 
peating  this  for  shots  at  different  stress  levels,  we  have 
also  been  able  to  obtain  the  pre-stress  dependence  of  the 
terminal  velocity.  The  results  for  a  glass  plate  are  shown 
in  Figure  29.  First  of  all#  the  velccities  indicated  are 
terminal  velocities.  We  could  not  determine  the  acceleration 
phase  of  the  crack  growth  since  our  measurements  could  not 
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begin  before  2  cm  radial  ' ;  tance.  Secondly,  there  were 
variations  in  velocities  between  c  =  0.7  and  1.7  km/sec, 
when  short  as  well  as  long  cracks  were  measured .  This  range 
is  indicated  in  Fig.  29. 

The  highest  crack  velocity  in  unstressed  glass  plate  is 
c  =  1.70  ±  .03  km/sec.  The  ratios  of  this  velocity  to  plate- 
compressional  and  :,hear  velocities  are  c/V  =  0,32  and 

r 

c/Vg  =  0.51  respectively.  For  the  glass  plate  stressed  to 

114  bars  tension,  the  dynamic  crack  velocity  was  c  =  1.34  ±  0.02 

km/sec  and  the  ratios  were  c/V  =  0.35  and  c/V  =  0.55.  The 

P  s 

intermediate  (45  bar)  tension  pre-stress  case  gave  values  the 
same  as  the  unstressed  case  within  the  experimental  error. 

It  would  be  useful  to  compare  our  values  of  crack  velo¬ 
cities  with  theoretical  prediction  as  well  as  other  experi¬ 
mental  results.  There  have  been  extensive  theoretical  studies 
of  crack  propagation.  Most  recent  reviews  can  be  found  in 
Liebowitz  (Vol.  I  and  II,  1968),  Tctelman  and  McEvillv  (1967) 
and  Cottrell  (1970).  Very  briefly,  a  necessary  condition  for 
the  propagation  of  an  elliptical  elastic  tensile  crack  is 
that  the  maximum  tensile  stress  at  its  tip  roach  the  theore¬ 
tical  cohesive  strength.  For  completely  brittle  and  elastic 
solids  where  the  elastic  limit  is  much  greater  than  the  cohe¬ 
sive  stress,  the  Griffith  relationship  can  be  applied.  Con¬ 
sidering  the  crack  to  be  a  thin  ellipsoid  with  semi-major  and 
minor  axes  1  and  h,  the  incremental  release  of  the  stored 
clastic  strain  energy  (Wp)  is  equal  to  l hr  i  nor-.-;  van l  al  increase 
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of  surface  energy  (Wg)  as  new  surface  is  created.  At  the  tip 
of  the  crack  there  is  a  strong  stress  concentration.  The 
transfer  of  this  elastic  energy  to  moving  the  adjacent 
points  apart  at  the  crack  tip  limi.ts  the  velocity  of  the  crack 
growth.  With  these  physical  models,  an  expression  can  be  de¬ 
rived  for  the  propagation  velocity  of  the  crack  (Tetelman  and 
McEvilly,  1967): 


=  0.38  V_  [  1  -  ~  ]l/i  «  0.38.  V_  [  1  "  8^1+V>YS  ^ 
P  Wt?  P  I  i 

E  TtZ*l 


where  c  =  crack  velocity,  V  =  compressional  velocity  (plate 

r 

velocity  in  our  case)  and  W  =  surface  energy  of  crack. 

Wg  =  strain  energy  in  the  crack,  ys  “  unit  surface  energy  and 
S  -  theoretical  strength,  l  -  half  crack  length  and  v  =  Poisson's 

ratio.  For  soda-lime  glasses,  for  example,  approximate  values 

3  2  '9  2 

are  yg  =  10  ergs/cm  m  and  Z  ~  10  dynes/cm  (Shand,  1961). 

From  these  it  is  clear  that  the  second  term  in  the  brackets 

is  fairly  small  and  the  maximum  theoretical  crack  velocity  is 

c  =  0.38  V  .  This  theoretical  value  is  in  agreement  with  some 

other  predictions  for  tensile  fractures  in  brittle  material 

where  a  value  of  about:  c  =  0.63  V  is  obtained  (Yoffe,  1951; 

w 

Mansinha,  1964).  There  have  been  other  theoretical  derivations 
relating  the  terminal  crack  velocity  to  Rayleigh  velocity 
(Baker,  1962;  Broberg,  1960;  Craggs,  1960;  Kachanov,  1961; 

Stroh,  1957).  In  these  cases  where  crack  velocity  coincides 
with  Rayleigh  velocity,  the  theoretical  values  differ  from 
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others,  unless  the  Rayleigh  velocity  is  taken  to  be  lover  at 
the  crack  tip. 

The  effect  of  pre-stress  in  the  material  is  to  increase 
stress  concentration  and  the  strain  energy  WE,  while  reducing 
the  bonding  energy  W  .  Thus  as  the  pre-stress  level  increases, 
WS/WE  approaches  zero,  and  the  crack  velocity  reaches  its 

asymptotic  value  of  c  +  0.38  V  . 

P 

The  agreement  between  our  measured  values  of  crack  velo¬ 
city  and  the  theoretical  values  of  eq.  16  is  very  good,  both 
in  absolute  value  and  pre-stress  dependence.  Ihe  observed 
values  are  slightly  lower  than  the  theoretical  values  based 
or.  perfectly  brittle  behavior.  This  is  expected  since  there 
is  always  some  plastic  deformation  in  the  fracture  process, 
and  soft  glasses  always  fail  to  reach  the  predicted  terminal 
velocities. 

There  have  been  other  measurements  of  crack  velocities 
in  glass  (Schardin,  1959} ,  and  those  are  in  excellent  agree¬ 
ment  with  ours.  The  crack  velocities  have  also  been  measured 
in  rock  by  a  high  speed  photographic  technique  similar  to 
ours  (Bieniawski,  1966).  In  gabbro  plates,  Bieniawski's  value 
for  terminal  fracture  velocity  was  c  =  1875  m/soc.  In  terms 
of  plate  velocity  V  ,  this  value  is  equivalent,  to  c  ~  0.31  V 

v  r 

in  excellent  agreement  with  our  value  under  zero  stress. 

The  implications  of  the  explosion -induced  crack  growth 
in  terms  of  the  radiation  pattern  of  seismic  waves  from  under¬ 
ground  explosions  ore  significant.  First  of  all,  •  ho  cracking 
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c nds  t.Vo.  re. g i  on  ,  and  a1  so  provides  fcho  means  of 

converting  some  explosion  energy  into  seismic  waves  of  compli¬ 
cated  radiation  pattern.  This  radiation,  however,  is  iso¬ 
tropic  in  tho  case  of  unstressed  media.  In  the  presence  of 
a  stress  field  or  tectonic  lineations,  the  cracks  grow  in  a 
preferred  direction  and  to  greater  lengths,  resulting  in 
coherent  radiation  of  tectonic  strain  energy  as  well  as  some 
converted  explosion  energy.  Even  under  ideal  experimental 
conditions,  however,  some  complexities  are  introduced  into 
the  radiation  patterns  due  to  the  presence  of  several  cracks, 
and  irregularities  in  growth.  These  may  explain  the  devia¬ 
tions  of  some  data  points  from  the  idealized  composite  radia¬ 
tion  patterns,  in  the  case  of  field  observations  such  as 
those  of  Figs.  4  and  8. 


IV.  DISCUSSION  AND  CONCLUSIONS 

So  far  in  this  paper  we  have  described  separately: 

(1)  the  source  properties  of  several  underground  explosions 
as  derived  from  the  analysis  of  seismograms ,  and  (2)  the  re¬ 
sults  of  a  series  of  laboratory  experiments  dealing  with  ex¬ 
plosive  sources  in  prestressed  media.  In  this  section  we  will 
combine  the  results  of  both  stiidies  to  explain  the  mechanisms 
of  Love  wave  generation  and  tectonic  strain  release  by 
explosions . 

The  main  results  of  field  data  analysis  and  the  associated 
laboratory  experiments  are  the  following: 

1.  All  explosions  detonated  in  hard  media  (granite, 
rhyolite,  tuff)  generate  some  Love  waves  at  or  in  the  immediate 
vicinity  of  the  source  (see  Table  1) .  Laboratory  experiments 
indicate  that  the  transverse  waves  are  generated  primarily 
because  of  the  existence  of  a  stress  field.  Even  in  brittle 
materials  such  as  glass,  under  unstressed  conditions,  few 
transverse  waves  are  generated  by  an  explosive  source  compared 
to  the  prestressed  case. 

2.  The  Love  and  Rayleigh  wave  radiation  patterns  can  be 
explained  by  a  composite  source  consisting  of  an  isotropic 
explosive  source  plus  a  double-couple  source.  It  was  assumed 


that  the  double-couple  component  was  horizontal  (pure  strike 
slip,  vertical  fault  plane).  With  the  available  scattered 
data,  it  is  not  possible  to  determine  uniquely  and  reliably 
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more  definite  source  models  specifying  all  tho  parameters 
(strike,  dip  angle,  slip  direction,  HUi«jnUmtn)  of  t1»r> 
function. 

3.  The  strength  of  the  double-couple  source  relative 

tc  the  explosion  (F  value)  increases  with  increasing  strength 
of  the  medium  as  shown  in  Table  1.  It  should  be  remembered 
that  (a)  F  is  period  dependent,  and  che  value  listed  in  the 
Table  is  applicable  at  T  =  15-20  seconds,  and  (b)  in  general 
the  F  value  listed  is  a  lower  bound  since  we  assumed  a  verti¬ 
cal  strike  slip  source  function  for  Love  wave  generation. 

This  is  the  most  efficient  radiator  of  Love  waves,  relative 
to  Rayleigh  waves.  Any  other  double-couple  source  model 
that  may  fit  the  observed  data  would  have  a  larger  F  value. 

4.  The  orientation  of  the  strike  of  the  double  couple 
in  general  is  in  agreement  with  the  strike  of  the  tectonic 
fiber  as  determined  from  the  existing  faults,  strikes  of  earth¬ 
quakes  and  the  aftershocks  following  the  large  explosions  such 
as  Benham  (Hamilton  and  Healy,  1969;  Kehrer,  1969;  McKeown 

and  Dickey,  1969;  Toksdz,  ejt  al .  ,  1965).  This  phenomenon  is 
in  agreement  with  what  was  observed  in  the  experiments  with 
an  explosive  point  source  in  prestressed  plates,  where  growth 
of  the  cracks  and  the  radiation  pattern  of  S-wavos  were 
.  governed  by  the  principal  stresses. 

5.  The  corrected  spectra  and  the  resulting  time  functions 
for  the  Love  and  Rayleigh  waves  showed  that  the  Love  wave 
spectra  were  relatively  richer  at  longer  periods  than  Rayleigh 
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wave  spectra.  In  fact/  Bilby  data  as  well  as  others  indicate 
that  the  spectra  of  Love  waves  generated  by  the  explosions 
are  similar  to  those  of  comparable  magnitude  earthquakes. 

Explosion  generated  Rayleigh  xaves,  on  the  other  hand,  have 
spectra  shifted  to  shorter  periods  (Toks&z,  et  al , ,  1965; 

Molnar,  et  al.(  1969;  Toksdz,  e±  al . ,  1970).  Our  one¬ 
dimensional  laboratory  experiments  also  show  that  the  time 
function  of  stress  relaxation,  generated  in  a  pre-stressed  rod,  is 
of  much  lower  frequency  compared  to  those  of  the  explosions  (see 
Figure  16) . 

6.  The  energy  of  the  double-couple  component  of  the 

source  holds  the  main  clue  to  the  source  of  this  energy  and 

the  mechanism  of  its  release.  Although  it  is  difficult  to 

compute  the  total  energy,  at  all  frequencies,  the  power  of 

the  double-couple  source  relative  to  the  explosion  at  T  -  20 

2 

seconds  is  proportional  to  F  .(Table  1) .  The  theoretical 
estimations  of  the  strain  energy  release  due  to  the  introduc¬ 
tion  of  a  cavity  in  the  pre-stressed  medium  are  hampered  by 
the  lack  of  direct  knowledge  of  pre-stress  (Press  and  Archambeau, 
1962;  Toksfiz,  et  ctl.  ,  1965;  Archambeau  and  S^mmis,  1970;  Smith, 
et  al. ,  1969).  The  laboratory  experiments  provide  some  answers. 

’  From  the  strained  rod  experiments  we  see  that  the  S-wavc  strain 
amplitude  increases  linearly  with  pvo-r'.  ress  (Figure  17)  as  ex¬ 
pected  theoretically.  The  calculations  carried  out  using  data 
from  unstressed  plate  experiments  shew  that  relaxation  around 
the-  cavity  alone  cannot  account  for  the  radiated  strain  energy. 
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Thin  combined  with  the  observed  G-wavc  radiation  patterns, 
requiring  a  growing  crack,  suggest  that  relaxation  associated 
with  extended  cracking  contributes  significantly  to  the 
radiation  of  strain  energy  in  pre-s tressed  media.  In  the 
case  of  field  testing,  relaxation  of  the  medium  occurs  both 
along  the  newly- formed  extended  cracks  (faults)  and  ah 
movement  along  pre-existing  fault  zones.  Strain  measure¬ 
ments  around  the  Nevada  Test  Site  show  this  (Dickey,  1969). 

8.  Based  on  field  studies  as  well  as  our  laboratory 
experiments  we  can  state  that  explosions  in  pre-stressed 
media  release  some  of  the  strain  energy,  and  the  mechanism 
of  the  release  of  this  energy  is  similar  to  that  of  an  earth¬ 
quake.  The  amount  of  tectonic  strain  energy  released  depends 
strongly  on  the  pre-existing  stress  level.  Explosion  induced 
cracking  controls  the  radiation  pattern  and  the  extent  of 
tectonic  strain  energy  release.  The  explosion  cavity  alone 
cannot  account  for  all  the  energy  release.  The  radiation 
patterns  of  seismic  waves  are  frequently  complicated  by  the 
growth  of  cracks  with  finite  velocities,  and  relaxations 
along  different  faults  and  cracks. 
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FIGURE  CAPTIONS 


Fig.  1. 


Fig.  2. 


Fig.  3. 

Fig.  4. 


Fig.  5. 


Fig.  6. 

Fig.  7. 


Long  period  seismograms  of  the  Bilhy  explosion  at  two 
LRSM  stations:  CP-CL  (Campo,  California,  A=482  km)  and 
FR-MA  (Forsyth,  Montana,  A=-1282  km).  Time  marks  are 
10  seconds  apart.  Z  indicates  vertical,  R  and  T  the 
horizontal  components,  respectively.  Nu:nbers  are 
relative  magnification  of  each  component.  Note  the 
absence  of  Love  waves  at  FR-MA. 

Love  and  Fayleigh  waves  from  the  explosion  Greeley, 
recorded  at  Kansas  City,  Missouri.  Z,  H  are  vertical 
and  horizontal  components,  respectively. 

Distribution  of  LRSM  stations  which  recorded  waves  from 
Biiby  in  North  America.  "Star"  indicates  shot  location. 
The  radiation  pattern  of  normalized  Rayleigh  waves  for 
Biiby.  The  points  are  the  ratio  of  the  amplitudes  of 
explosion-generated  to  collapse-generated  Rayleigh 
waves.  The  curve  is  the  theoretical  radiation  pattern 
for  a  composite,  source  consisting  of  an  explosion  and 
an  orthogonal  double-couple. 

Rayleigh  wave  amplitude  response  of  the  layered  medium 
to  an  explosive  source  near  the  surface,  with  impulsive 
time  function. 

Love  wave  amplitude  response  of  the  medium  to  a  double¬ 
couple  source  near  the  surface. 

Love  wave  to  Rayleigh  wave  vertical  component  Fourier 
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Fig,  3.  Arnpl  i  f.udo  ratios  of  the  explosion-generated  Love  and 
Rayleigh  waves  (Ug/URz)  as  a  function  of  azimuth. 

The  theoretical  curve  is  for  the  composite  source 
described  in  Fig.  4. 

Fig.  9.  Rayleigh  v/ave  pulses  at  tv/o  stations  plotted  from 
digitized  data. 

Fig.  10.  Fourier  amplitude  spectra  of  the  Rayleigh  v/ave  pulses 
shewn  in  Fig.  9. 

Fig.  11.  Ground  displacement  spectra  obtained  from  those  of 
Fig.  10  after  correction  for  instrument  response. 

Fig.  12.  Amplitude  spectra  of  source  time  function  after  cor¬ 
rection  for  propagation  effects.  Circles  indicate 
data  uncorrected  for  attenuation?  triangles  ind’eate 
the  data  corrected  taking  Q  =  100. 

Fig.  13.  The  Rayleigh  v/ave  source  time  function  for  Bilby  at 
the  boundary  of  the  seismic  zone  where  strains  are 
infinitesimal. 

Fig.  14.  Amplitude  spectra  of  source  time  function  for  Love 
waves  after  correction  for  propagation  effects. 
Circles  are  uncorrected;  triangles  are  corrected  for 
attenuation  taking  Q  =  100. 

Fig.  15.  Laboratory  set-up  for  one-dimensional  experiment. 

Stress  waves  from  an  explosive  source  in  a  rod  under 
torsion  are  recorded  by  strain  gauges. 

Fig.  16.  Strain  gauge  seismograms  from  one-dimensional  experi¬ 
ment  at  different  pre-stress  values. 
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Fig.  17.  .'.aptitudes  of  P  and  S  vmvos  as  a  function  of  pre-strain 

Fig.  18.  Experimental  set-up  for  the  two-dimensional  dynamic 
photoelastic  studies  of  explosive  sources  in  pre¬ 
stressed  media. 

Fig.  19.  Isochromatics  and  isoclinics  from  a  shot  in  unstressed 
1  "  . 

-j-g—  thick  plexiglas.  Interframe  time  4.16  yscc. 

Fig.  20.  Theoretical  isochromatics  and  isoclinics  for  P  and  S 
waves . 

Fig.  21.  Theoretical  pulse  shapes  as  a  function  of  distance  for 
radial  stress  and  displacement  and  maximum  shear  stress 
from  a  cylindrical  source  in  a  plate,  r  is  distance 
in  units  of  source  radius. 

1" 

Fig.  22.  Isochromatics  and  isocimics  from  a  shot  in  a  ^  thick 
glass  plate  stressed  under  tension  (45  bars).  Inter¬ 
frame  time  2.08  pscc.  Note  the  growth  of  cracking. 

Fig.  23.  Dynamic  radiation  patterns  of  P  and  S  waves  from  ex¬ 
plosive  sources  in  pre-stresscd  plexiglas  and  glass 
plates  under  tension,  as  measured  by  strain  gauges. 

Fig.  24.  Maximum  strain  energy  ratios  of  P  and  S  waves  versus 
azimuth  from  tensile  axis. 

Fig.  25.  Theoretical  P  and  S  wave  radiation  patterns  from  a 
tensile  fault  plus  explosion  in  glass.  Each  curve 
individually  normalized  to  its  extreme  value.  «,  8 
are  compressional  and  shear  velocities  in  the  plate, 
c  =  crack  propagation  velocity.  the  ratio 

of  the  strength  of  the  isotropic  component  to  tensile 
crack  component. 
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TABLE  2 

Observed  S-wave  and  calculated  maximum  available  energies 
for  explosions  in  pre-stressed  plexiglas  and  glass  plates. 


Maximum  available 


Shot  No. 

Medium 

Tensile 

pre-stress 

(bars) 

Observed  S- 
v/ave  energy 

strain 

a  = 

energy 

1  cm 
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Plexiglas 

163 

49.3 

5 

x  10 

2.84 

,n5 

X  JlO 
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Plexiglas 

163  (?) 

3.72 

x  105 

CO 

• 
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III.  PILE  DRIVER  AND  GREELEY:  TECTONIC  UTILITY  AND  DANGERS 
OF  NUCLEAR  EXPLOSIONS 


Published  in  Science,  173 , • 230-233 ,  July  16,  1971. 


-  123 


ABSTRACT 

The  tectonic  strain  energy  released  by  several  under¬ 
ground  nuclear  explosions  has  been  calculated  through  an 
analysis  of  seismic  surface  waves.  The  proportionally 
great  amount  of  energy  released  in  the  case  of  certain  events 
suggests  both  the  possible  uses  for  as  well  as  hazards  of 
underground  testing. 
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Seismic  waves  generated  by  underground  nuclear  explosions 
quite  often  indicate  complications  arising  from  the  interaction 
of  the  explosive  source  with  the  surrounding  inhomogeneous 
medium.  The  clearest  evidence  of  such  an  interaction  is  the 
generation  of  horizontally  polarized  SH  and  Love  waves  by  many 
explosions.  The  radiation  patterns  of  these  waves  and  the 
mechanisms  of  their  generation  have  been  a  subject  of  a  number 
of  studies  (1) .  It  is  generally  agreed  that  an  explosion  deto¬ 
nated  in  a  pre-stressed  medium  radiates  some  ambient  strain 
energy  because  of  relaxation  around  the  cavity  and  extended 
cracks.  The  extent  of  this  relaxation  and  the  amount  of  energy 
released,  however,  have  not  been  accurately  determined. 

The  radius  of  the  zone  of  relaxation  and  the  amount  of 
tectonic  strain  energy  release  are  two  characteristics  of 
nuclear  explosions  which  are  of  extreme  importance.  The  first 
consideration  which  depends  upon  a  knowledge  of  these  charac¬ 
teristics  is  the  feasibility  of  utilizing  underground  explosions 
for  earthquake  control.  Explosions  might  be  effective  in 
periodically  releasing  accumulating  tectonic  strain  energy 
in  active  areas,  thus  possibly  preventing  major  earthquakes. 

The  second  consideration  concerns  the  safety  of  large  under¬ 
ground  explosions.  An  explosion  which  releases  very  large 
amounts  of  strain  energy  would  have  the  same  effect  as  that 
of  a  large  earthquake.  If  larger  devices  are  tested  in  new 
locations  where  ambient  stress  patterns  are  not  known,  the 
strain  release  problem  becomes  more  critical.  As  a  third 
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point,  the  release  of  large  amounts  of  pre-existing  strain 
energy  may  affect  some  (but  not  all)  of  the  seismic  discrimi¬ 
nation  methods  between  earthquakes  and  explosions. 

The  cases  that  have  been  studied  to  date  indicate  that 
some  strain  energy  was  released  by  most  nuclear  explosions 
detonated  in  relatively  hard  media  (2) .  However,  the  strain 
energy  released  was  less  than  the  equivalent  seismic  energy 
of  the  explosions  themselves  in  all  but  one  case.  The  ex¬ 
ception  was  Hardhat,  detonated  in  granite.  The  excessive 
strain  energy  release  in  this  case  was  attributed  to  the 
possible  triggering  of  an  earthquake  (3) . 

In  this  paper  we  will  report  on  two  other  Nevada  Test 
Site  explosions  (Pile  Driver  and  Greeley)  which  generated 
large  Love  waves  and  released  significant  amounts  of  tectonic 
strain  energy.  Data  from  these  events  was  obtained  from 
LRSM  (Long  Range  Seismic  Measurements) ,  WWSS  (World  Wide 
Standard  Seismograph) ,  and  CSS  (Canadian  Standard  Seismograph) 
system  stations  located  throughout  North  America.  The  long 
period  records  were  digitized  and  filtered.  The  two  hori¬ 
zontal  components  at  each  station  were  rotated  to  radial  and 
tangential  directions  with  respect  to  the  explosions  in  order 
to  separate  the  Love  waves.  The  resulting  tangential  compo¬ 
nent  (Love  waves)  and  the  vertical  component  (Rayleigh  waves) 
were  then  Fourier  analyzed  to  obtain  the  amplitude  spectra. 

In  Figure  1,  long  period  seismograms  of  three  explosions 
are  illustrated.  Tan,  detonated  near  the  other  events,  has 
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been  included  as  an  example  of  a  '  y:>ical  explosion  which  did 
not  generate  extensive  Love  wave-, .  in  this  study  both  the 
radiation  patterns  of  Raylov-h  v  •»<-;•  and  Love  to  Rayleigh 
wave  amplitude  ratios  were  ux-c-J  !  a  I'oterrr.jne  the  strain 
energy  release. 

An  explosion  with  associated  tectonic  strain  release  can 
be  represented  by  a  composite  source  consisting  of  an  isotropic 
explosion  component  and  a  double-couple  strain  release  compo¬ 
nent  (4) .  .  Assuming  that  the  time  delay  between  the  explosion 
and  the  strain  release  as  well  as  the  differences  between 
Love  and  Rayleigh  wave  time  functions  are  negligible,  expres¬ 
sions  for  Rayleigh  and  Love  wave  far  field  displacements  from 
such  a  composite  source  are  given  by 


Ur(u)  =  CR(w) Cl  +  F  sin  20]  exp[i (et-kRr) ] 
cos  20  exp[i  (ut-k^r)  ] 


(1) 


where  CR  and  are  functions  of  frequency,  distance,  and  the 
medium,  w  is  frequency,  t  is  time,  kR  and  kL  are  the  Rayleigh 
and  Love  wave  numbers,  and  r  is  the  radial  distance  from  the 
source.  0  is  the  azimuthal  orientation  of  the  double-couple 
and  F,  the  parameter  of  greatest  interest,  is  the  strength  of 
the  double-couple  relative  to  the  explosion. 

If  the  differences  between  Love  and  Rayleigh  wave  attenua¬ 
tion  and  source  time  functions  arc  considered  negligible,  the 
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following  expression  is  obtained  for  the  Love  over  P.ayleigh 
wave  amplitude  ratio 


Ki 

i"3 


Clr(w)  F  cos  20 
1  +  F  sin  26 


(2) 


where  C^R  is  a  function  of  frequency  and  the  medium.  An  auto¬ 
matic  error  scheme  was  applied  to  fit  the  observed  data  to 
this  formulation  (5) .  A  grid  was  formed  between  the  two  para¬ 
meters  F  and  9,  and  the  error  between  the  data  and  the  theore¬ 
tical  combinations  of  F  and  0  was  contoured. 

Figure  2  shows  the  best  fitting  Love  over  Rayleigh  wave 
radiation  patterns  for  Pile  Driver  and  Greeley  along  with  the 
observed  data.  The  parameter  F  determined  for  these  two  events 
was  3.2  and  1.6.  respectively.  In  Figure  3a,  the  best  fitting 
Rayleigh  wave  pattern  for  Pile  Driver  is  shown  with  the  observed 
data  and  some  actual  seismograms.  In  addition  to  amplitude, 
phase  (polarity)  has  been  included.  The  polarity  reversals 
predicted  by  our  best  fitting  theoretical  model  (since  F  >>  1) 
are  indeed  observed  when  the  Pile  Driver  data  is  compared  to 
that  of  Tan. 

The  Rayleigh  wave  radiation  pattern  for  Greeley  is  shown 
in  Figure  3b.  In  this  case,  the  strength  of  the  strain  release 
component  was  close  to  that  of  the  explosion,  and  the  amplitudes 
of  the  inverted  lobes  were  therefore  small.  In  spite  of  this, 
a  very  clear  phase  reversal  can  be  observed  relative  to  the 
Boxcar  event  as  seen  at  LASA  from  the  data  of  Filson  (6) . 
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The  source  model  for  the  Hard hat  event  was  studied,  by  a  simi¬ 
lar  technique  and  described  in  a  previous  publication  (7). 

The  results  were  almost  identical  to  Pile  Driver,  and  clear 
polarity  reversals  could  be  observed. 

Since  the  radiation,  patterns  of  P  waves  arc  similar  to 
those  of  Rayleigh  waves,  we  investigated  these  body  waves  to 
determine  any  polarity  reversals.  The  short  period  records, 
however,  appeared  to  be  dominated  by  the  explosive  component 
end  near-source  irregularities.  Only  on  the  long  period  com¬ 
ponents  of  3  stations  did  P  waves  from  Greeley  appear  reversed 
relative  to  Boxcar.  In  the  case  of  Pile  Driver,  the  long 
period  P  waves  were  far  too  small  to  be  useful. 

The  source  data  for  the  three  explosions  of  concern, 
Hardhat,  Pile  Driver,  and  Greeley,  as  well  as  that  for  some 
other  events,  are  summarized  in  Table  1.  In  addition  to 
double-couple  strengths,  we  have  computed  the  ratio  of  surface 
wave,  energy  of  the  double-couple  component,  Etcct  >  to  the 
energy  of  the  explosion.  This  ratio  (8)  is  closely  approxi¬ 
mated  by 

Etect  /  Ee*p  *  V3  F2. 

Several  conclusions  can  be  derived  from  these  data. 
Generally  (in  the  Nevada  Test  Site) ,  explosions  in  harder 
media  such  as  granite  release  more  tectonic  strain  energy 
than  explosions  in  softer  formations  such  as  tuff  or  alluvium. 
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These  latter  media  have  low  rigidities  and  cannot  accumulate 
shear  strain  energy  because  of  plastic  deformation  under  stress. 
Salt  is  a  very  good  example  of  such  a  plastic:  medium.  None 
of  the  explosions  in  salt  generated  measurable  Love  waves. 

Not  all  explosions  in  a  given  medium,  however,  have  propor¬ 
tionally  the  same  strain  release  associated  with  them,  as 
indicated  by  Shoal  versus  Pile  Driver  and  Hardhat.  Thus  in 
addition  to  a  rigid  medium,  high  ambient  tectonic  stresses 
are  important  in  determining  the  strain  energy  release  by  an 
explosion.  The  similarity  of  the  Pile  Driver  and  Hardhat 
results  indicate  that,  1)  the  regional  tectonic  framework 
controls  the  radiation  pattern,  2)  in  a  given  medium  and 
area,  each  explosion  releases  strain  energy  proportional  to 
its  own  explosive  energy,  3)  the  spontaneous  energy  release 
comes  from  the  vicinity  of  the  shot  rather  than  from  an  earth¬ 
quake  triggered  at  a  distance  of  ten  kilometers  or  more. 

Strain  adjustments  at  such  distances  appear  to  take  place 
gradually  as  indicated  by  the  spread  of  aftershock  activity 
that  follows  some  large  explosions  (9). 

The  implications  of  these  results  in  terms  of  the  utility 
of  nuclear  explosions  for  the  release  of  tectonic  stresses 
and  earthquake  control  are  encouraging-,  Since  an  explosion 
can  release  accumulated  strain  energy  considerably  greater 
than  its  own  (more  than  ten  times  greater  in  the  cases  of 
Pile  Driver  and  Hardhat)  in  the  vicinity  of  the  shot  point, 
such  events  could  be  used  to  relax  small  regions.  One 


Sr 


possible  application  would  bo  to  do-activato  a  site-  in  a  seis- 
micclly  active  region  prior  to  the  construction  oi  some  project 
such  as  a  nuclear  power  plant,  however,  at  the  present  we  do 
not  know  how  long  it  takes  for  the.  tectonic  stress  field  to 
build  up  to  its  initial  level  following  an  explosion.  The 
Hardhat  and  Pile  Driver  shot  locations  were  only  about  a  half 
kilometer  apart.  Yet  Pile  Driver,  detonated  four  years  later, 
radiated  proportionately  the  same  amount  of  strain  energy. 

This  does  not  necessarily  mean  that  strain  energy  reaccumulatec 
at  the  Hardhat  site  in  only  four  years.  Since  the  yield  of 
Pile  Driver  was  10  times  greater  than  that  of  Hardhat,  the 
stress  relaxation  (primarily  due  to  cracking)  most  likely  ex¬ 
tended  beyond  the  relaxation  zone  cf  the  h'ardhat  explosion. 

Thus  tectonic  strain  energy  was  released  from  a  larger  volume 
in  the  case  of  Pile  Driver.  The  use  of  explosions  to  release 
tectonic  strain  energy  in  areas  which  are  already  developed, 
however,  is  clearly  not  feasible. 

For  the  same  reasons  that  explosions  could  prove  useful/ 
nuclear  testing  in  regions  of  high  ambient  stresses  might 
have  serious  consequences.  Since  the  medium  properties  and 
not  the  shot  yield  control  the  proportion  of  energy  release, 
a  large  explosion  could  release  large  amounts  of  tectonic 
strain  energy.  For  example,  with  the  energy  factor  of  Pile 
Driver  or  Hardhat,  a  one-megaton  explosion  could  release  energy 
equivalent  to  that  of  a  magnitude  Mg  =  6.3  earthquake.  For 
a  ten  megaton  explosion  this  could  be  equivalent  to  a  magnitude 
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Mg  =  7.2  earthquake,  indeed  a  potentially  destructive  earth¬ 
quake.  Thus  in  any  testing  program,  those  factors  must  be 


taken  into  account,  and  hard  media  should  be  avoided  unless 


the  ambient  stress  levels  are  known  to  be  low. 
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figup.l  c tj  tiokc 

Fig.  1.  Long  period  filtered  seismograms  illustrating  relative 
generation  of  Love  waves  by  three  N.T.S.  explosions  recorded 
at  Weston,  Massachusetts.  Love  waves  appear  on  the  tangential 
components  (T)  about  two  minutes  earlier  than  the  Rayleigh 
waves  on  the  vertical  components  (Z).  Note  the  large  amplitude 
of  Love  waves  relative  to  Rayleigh  waves  for  Pile  Driver  and 
Greeley  explosions. 

Fig.  2.  Love  over  Rayleigh  wave  amplitude  radiation  patterns 
for  Pile  Driver  and  Greeley.  For  theoretical  curves  (solid 
lines),  the  following  parameters  were  used:  orientation  (0) 
and  strength  of  double  couple  (F) ,  340°  from  north  and  3.2 
respectively  for  Pile  Driver,  and  355°  and  1.6  for  Greeley. 

Data  symbols:  B  WWSS  stations;  A  LRSM  stations;  $  CSS 
stations . 

Fig.  3.  Rayleigh  wave  radiation  patterns  for  explosive  com¬ 
ponent  (circles),  strain  release  component  (broken  curves), 
and  composite  source  representing  the  theoretical  models  (solid 
curves,  pax*ameters  those  of  Fig.  2).  The  polarity  (phase)  of 
Rayleigh  waves  in  different  quadrants  of  the  radiation  pattern 
is  indicated  by  (+)  and  (-).  (a)  Pile  Driver  with  seismograms 

compared  to  those  of  Tan.  Pile  Driver  traces  are  below  those 
of  Tan.  Dashed  line  traces  have  their  polarities  reversed. 

Note  the  perfect  match  of  wave  shapes,  with  polarities  reversed 
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at  JP-AT  and  B02  us  predict.*..:1  by  the  model.  (b)  Greeley  \:i 
seismograms  compared  to  those  of  Jin  I  f  bonk  (COL  and  OXF)  and 
Boxcar  (LASA) .  The  lower  tracer  are  from  Greeley.  Polarit 
reversal  is  very  clear  at  LASA.  Data  symbols:  A  LP.S.M 
stations;  A  LRSM  stations  where  reversed  polarity  v/ns 
observed;  S  WWSS  stations;  3  CSS  stations. 
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